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Abstract 

We prove van der Waals - London's law of decay of the van der Waals force for a 
collection of neutral atoms at large separation. 

1 Introduction 

Van der Waals force between atoms and molecules plays a fundamental role in quantum 
chemistry, physics and material sciences. For instance, it defines the chemical character of 
many organic compounds and enable geckos - which can hang on a glass surface using only 
one toe - to climb on sheer surfaces (see [Wiki j ) . It explains why water condenses from vapor 
as well as many properties of molecular compounds, including crystal structures (e. g. the 
shapes of snowflakes), melting points, boiling points, heats of fusion and vaporization, surface 
tension, and densities. It forces gigantic molecules like enzymes, proteins, and DNA into the 
shapes required for biological activity (see [Sen] ). 

A microscopic explanation of this force was given by F. London soon after the discovery 
of quantum mechanics and was one of its early triumphs. This heuristic explanation showed 
that this force has a universal behavior at large distances - it decays as the inverse sixth power 
of the distance between atoms. This behavior was confirmed in [LTj by proving - through 
a sophisticated test function construction - an upper bound. Our goal in this paper is to 
provide a complete proof of the van der Waals - London decay law for atoms. 

Let — e and m denote the electron charge and mass. Consider a system of M multielectron 
atoms which we call a molecule though we do not assume binding between atoms. In the 
units where h = 1 the hamiltonian of the system is 

1 . ^ e 2 ^ , ^ e 2 ^ 1 . .'A ZiZj* 



2m 1 ^ \Xi -yj\ ^ \Xi -Xj\ ^ 2mj yj \yi - y-\ 

Here N is the total number of electrons, Xi,yi € M 3 denote the coordinates of the electrons 
and the nuclei, respectively, eZj is the charge of the j-th nucleus and rrij is the mass of the 
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j-th. nucleus. We consider a system of neutral atoms so we must have X^=i %j = N. The 
operator H mo \ acts on the subspace of the space L 2 (M?( N+M )), which accounts for the 

fact that the electrons are identical particles and are fermions and therefore they obey the 
Fermi-Dirac statistics. Also, possibly, some of the nuclei are identical and obey either the 
Fermi-Dirac or Bose-Einstein statistics (more details are given below). For M = 1, H mo \ is 
the Hamiltonian of the atom with the nucleus of charge eZ = eZ\. 

Let E m)n denote the ground state energy of the ion with a nucleus of charge eZ m and 
Z m — n electrons, so that E m) Q = E m is the ground state energy of the atom corresponding 
to the m-th nucleus of charge Z m , and let y = {y\, ...,yM) be the collection of the nuclear 
co-ordinates. The interaction (or cohesive) energy W(y) between atoms in this system is 
defined as 

W(y):=E(y)-E(oo), (1.1) 

where E{y) is the ground state energy of the system with positions of the nuclei fixed, i.e. 
the ground state energy of the hamiltonian 

N 1 M p 2 7 . 1 ' JV p 2 M/ 7? z 

H N (y) = " E r^-r) + E r^H + E v^~v ^ 

2m ,\Xi — yi\ .\xi — Xj\ . \Vi — Vjl 

acting on the subspace %f crm i of the space L 2 (U. 3N ), which accounts for the Fermi-Dirac 
statistics of electrons, and E(oo) is the energy of the system with the atoms infinitely far 
from each other, 

M M 

E(oo) =mm{J2 E m,n m \ E n ™ = °}" ( L3 ) 

m=l m=l 

(Of course, one expects that E(oo) = z~2m=i ^m, see the discussion below.) As we will see 
below E(y) as well as E m are isolated eigenvalues of finite multiplicities. Also by the HVZ 
theorem below, W6 Ccin retain in (|1 .3p only terms for which E m ^ nrn are isolated eigenvalues. 

The hamiltonian (|1.2p is called the Born-Oppenheimer hamiltonian. It arises as the central 
object in the key technique in solving the eigenvalue problem for i/ mo i, which is called the 
Born-Oppenheimer approximation, and which plays an important role in quantum chemistry 
(for example minima of its ground state energy E{y) determine shapes of molecules). In this 
approximation the ground state energy E(y) (or the energy of an excited state) of H^(y) is 
considered as the potential energy of the nuclear motion, which leads to the Hamiltonian 

M 



One expects that due to the fact that the ratio of the electron and nuclear mass being very 
small, the eigenvalues of -ff nuc i give a good approximation to the eigenvalues of H mo \. For 
rigorous results on the Born-Oppenheimer approximation see original articles [CDS! IKMSW} 
iLeL] [Hagj [Hag2| IPST] and the textbook [O^] . 



One can define the interaction energy to any order in in the electron to nuclei mass ratio 
(see Appendix|B]), but this will produce only correspondingly small corrections to our results. 
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It is expected, after van der Waals, that W(y) is a sum of pair interactions, Wij, which 
are attractive and decay at infinity as — \yi — Vj\~ 6 - More precisely, one expects that 

1,M 4 1,M 4 

^-£^ + o ( £^), (1.4, 

where o~ij are positive constants depending only on the parameters of the pair of atoms i, j, 
provided that min{|yj — y~j\ : 1 < i < j < M} is large enough. 

We formulate a property of many body systems playing an important role below 

(E) EZi < E-ii E i>ni , V(m, . . . n M ) = 0. Ei M > 0- 

Property (E) depends on the symmetry type <r, defined below, since the ground state 
energies depend on the underlying spaces on which the hamiltonians are defined. Experiments 
and numerical computations (see below) show that Property (E) holds for all elements for 
which it was tested, however, theoretically, it proven only for a system of several hydrogen 
atoms. If it holds, then E(oo) = J2m=i Em- 

The simplest symmetry type for fermions of spin ^ is the one corresponding to totally 
antisymmetric functions, or, what is the same, to the greatest possible spin. We begin with 
formulating our restlts in this special case. 

Theorem 1.1 (van der Waals law; highest spin). Assume that the Hamiltonian H^{y) acts 
on the space Ha of purely antisymmetric functions and assume Condition (D) stated below. 
Then for large distances between atoms (jl.4p holds for some constants o~ij > depending on 
the nature of the atoms i,j, if and only if Property (E) stated above holds. 

The theorem above is a corollary of Theorem 11.21 stated below, dealing with general 
symmetry types. 

Now we define the physical state space, %f er mi) of the Born-Oppenheimer molecule. Since 
electrons are identical particles and are fermions of spin ^, the state space of the system of 
N electrons is the space 

N 

/\(L 2 (M 3 )®C 2 ) 
l 

of I?— functions, 9(xi, s±, . . . , x??, sjv), of co-ordinates, x\, . . . , Xn, and spins, Si,...,Sjv 
(with Sj € {— = l,...,iV) that are antisymmetric with respect to permutations of 
pairs (xi,Si). The space %f C rmi is the subspace of L 2 (M. 3N ) given by the projection of the 
space Af (L 2 (M 3 ) (8) C 2 ) onto the L 2 -funct ions of the co-ordinates alone, 

N 1 1 

"Hfcrmi := {(X, *>spin|* € /\(L 2 (R 3 ) C 2 ), X = {"g. ^ ~* C} ' 

where (%, *) sp i n := E si ,..., sjve {-I,i} x( s i> ■ • ■ > sjv)*(^i, H, ■■■,xn> s n ). 

We relate this space to irreducible representations, T a , of the group Sn, of permutations 
of N indices. Consider the unitary representation T : Sn —> U(L 2 (M. 3N )) (unitary operators 
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on I?) of S N on the space L 2 (R 3N ), in detail, n G S N -> T„ G C/(L 2 (M 3JV )), with the unitary 
operators given by 

(^^(xi,...,^) = *(a; n -i(i),a;^-i(2),...,a;^-i(j\r)). (1.5) 
Then the space "Hf erm i can be written as 

H ieimi = J2K a , (1-6) 

where a runs over irreducible representations of the group Sjv corresponding to at most 
two-column Young diagrams and "H 17 is the subspace of L 2 (M. 3N ) on which this representa- 
tion reduces to multiple of the irreducible representation of the type a (see Section 17.11 for 
definitions and details). We call irreducible representation labels a the symmetry types. 

The fact that the electrons are identical particles is expressed in the property that H^{y) 
commutes with the permutations 

H N (y)T„ = T^H N {y), Vvr G S N , (1.7) 

and therefore the subspaces T~L a are invariant under Hj^(y). This allows us to introduce the 
ground state energy of the system for the symmetry type a by 

E°(y) = Ma(H N (y)\ n *). 

To define E a (oo) we need some notation. Let a be the decomposition of the molecule into 
neutral atoms and/or ions, S(a) C SV be the subgroup of SV consisting of the permutations 
that keep the clusters of a invariant and let H a be the sum of the corresponding atomic or ionic 
Hamiltonians (see the next section for precise definitions). Fix an irreducible representation 
T u of the group SV. Then the restriction T a \g^ is not necessarily irreducible. Therefore, 
there exists a family, I a , of irreducible representations of S(a) such that 

where T° is the irreducible representation of S(a) of the type a. The representations a G I a 
are called induced representations and we write a -< a. (For hydrogen atoms the group S(a) 
is trivial and the construction of induced representations should be omitted.) We write 

a a if a -< a and inf a(H°) = mininf a(H?), 

b,0 b 

where H% denotes the restriction of H a onto the multiple of the irreducible representation of 
type a). We now denote the lowest energy of separated atoms and/or ions for the symmetry 
type a as 

£ ct (cxd) = min infa(^). 
The interaction energy for the symmetry type a is now defined as 

W a (y) := E a {y) - E%<x). 



Finally, we state the Condition (D) of Theorem 11.11 and of the theorems below: 
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(D) For each a -<~< a the ground state subspace of H" consists only of one copy of the 
irreducible representation of type a. 

One expects that for every symmetry type a, the ground state subspace consists of a single 
copy of the irreducible representation of the symmetry group, but proving this is an open 
problem. For a system of several hydrogen atoms Condition (D) follows from the fact that 
S(a) is the trivial group and from Perron- Frobenious argument (see below). 

The following theorem gives the van der Waals law for fixed symmetry types. 



Theorem 1.2 (van der Waals forces for a fixed symmetry type). Assume Condition (D) 
below. Then for every symnn 
( defined in (|A.4j) ) such that 



below. Then for every symmetry type a and every a -<-< a, there exist positive constants 



where 



W a (y) = min W a > a {y), 

1,M a,a l,M 

if and only if Property (E) reinterpreted for the symmetry type a holds. 

For a collection of hydrogen atoms Condition (D) and Property (E) are shown below to 
hold and therefore the van der Waals law is always valid for such a system. 

As was mentioned above, using an intricate test function Lieb-Thiring |LT| have proven 
an upper bound for the interaction energy with a better control of how large interatomic 
separations should be. These bound also works for forces between certain molecules. When 
a corresponds to a Young diagram with one column (completely symmetric representation) 
Theorem 11.21 gives Theorem 11.11 For a corresponding to a Young diagram of at most two 
columns the ground state energy of Hj^(y) on %f crm i is 

E(y) = mmE a (y). 

a 

Let do be a Yonge diagram of at most two columns for which E(y) = E a °(y). Then The- 
orem 11.21 for the specified o"o gives the interaction energy of the system. (Of course, if the 
interatomic distances are not very large it might happen that the energy surfaces for different 
symmetries cross and we have to take W(y) = E a ° (y) — E ai (oo) where o~\ ^ o~o is the diagram 
that minimizes E a (oo).) 

Theorem 11.21 describes the van der Waals force at a pairwise large separation between 
the atoms. For intermediate distances, the van der Waals - London law is modified due to 
overlapping between electron clouds of the atoms and for small distances, the van der Waals 
forces are repulsive (the energy is positive) as follows from the rough estimate 

X > M e 2 7 7 . 



for some constant C independent of y, implied by the bound ^ jf™ < —aA Xn + f3, valid 
for any a > and a corresponding (3 > 0. Often the interaction energy for two atoms 
(M = 2) is modeled by the Lennard-Jones potential Wlj(v) = i — rn — i — ^ — m or by 
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the Buckingham potential Wb(u) = e~ c \ yi ~ V2 \ — r^-j^p where the constants a,b,c,d are 
determined experimentally. 

The van der Waals-London law is expected to be true for the case of molecules that are 
not polarized (including the noble gases and methane). If the molecules have dipole moments, 
then one expects a third power law to be true. This can be proven by the techniques developed 
in this paper. 

For a system of several hydrogen atoms Property (E) is proven in 

Proposition 1.3 (Property (E) for a system of several hydrogen atoms). For a system of 
several hydrogen atoms Property (E) holds and consequently van der Waals law holds always. 

Proof. Let E^ denote the ground state energy of the hydrogen atom and let E^ be the 
ground state energy of the hydrogen ion with charge —me (where m > 1), i.e. the lowest 
eigenvalue the Hamiltonian 

m+I l,m 2 

j=l i<j 1 J 1 

2 

where Hj = — A Xj — j|-y is the Hamiltonian for the hydrogen atom in the j-th coordinate. 
Property (E) for a system of several hydrogen atoms is reduced to the property that for any 
m > 1, £W satisfies 

E^ > (m + l)EW. (1.9) 

Let := max{m' < m| H^ m ) has a ground state} and let tp^ m ^ be the ground state 
of £f( m *) corresponding to E^ m *\ (We know from |Hij that > 1 for m > m,.) By 
the definition of m* we have E^ = E^ and therefore E^ = (^ rn *\ H^^ m ^). Since 
Hj > E for any j, we have E^ > (m+l)E^+5, with 5 := (V> (m *\ El<7* j^-\i> {m * ] ) ■ □ 

Below we show that Property (E) follows from the following stronger condition 

(E') For any two nuclei i and j, i ^ j, in our system and for m + n < Zi, —n < Zj, mn > 
0, n ^ 0, we have that 

Note that Property (E) for a system of hydrogen atoms with electron statistics follows 
from Property (E) without statistics. Indeed, E^ remains the same whereas E^ increases 
if the statistics is taken into account. 

The meaning of Property (E') is that ionization energies of atoms are greater than the 
electron affinities, where, in a standard terminology, the n-th ionization energy (n > 1) of an 
atom is the energy required to remove an electron from its n — 1-ion and the n-th electron 
affinity is the energy required to remove an electron from its — n ion. The table below, taken 
from [MA], gives first ionization energies and first electron affinities of atoms. It shows that 
the first ionization energies are always much larger than the first electron affinities. 
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Atomic number 


Element 


first Ionization 
energy (kcal/mol) 


Electron affinity 
(kcal/mol) 


1 


H 


313.5 


17.3 


2 


He 


566.9 


_ 


3 


Li 


124.3 


(14) 


4 


Be 


214.9 




5 


B 


191.3 


(7) 


6 


C 


259.6 


29 


7 


N 


335.1 


_ 


8 





314.0 


34 


9 


F 


401.8 


79.5 


10 


Ne 


497.2 


_ 


11 


Na 


118.5 


(19) 


12 


Mg 


176.3 




13 


Al 


138.0 


(12) 


14 


Si 


187.9 


(32) 


15 


P 


254 


(17) 


16 


S 


238.9 


47 


17 


CI 


300.0 


83.4 


18 


Ar 


363.4 


(16) 


19 


K 


100.1 




20 


Ca 


140.9 




21 


Sc 


151.3 




22 


Ti 


158 




23 


V 


155 




24 


Cr 


156 




25 


Mn 


171.4 





(Values in parentheses are estimated by quantum-mechanical calculation and have not been 
verified experimentally.) 

It is experimentally verified that the higher (second third and so forth) ionization energies 
are bigger than the first and it is expected that the higher electron affinities affinities are lower 
or zero. (Only ions with at most two extra electrons - i.e. of the charges at most — 2e - are 
observed experimentally. That nuclei can bound only finite number of electrons was proven 
in [Rusl [Sig] (with the bound < 17 Z on the number of the extra electrons). It was shown in 
[LSST] , with some improvements by [SSSt IFSlj , that asymptotically ions are neutral and in 
[Lieb| that the maximal number of extra electrons < Z + 1. The latter bound was recently 
improved to < 0.22Z + 3Z 1 / 3 in [Phanj . Moreover, it was shown in [Solll [S0I2] that in the 
Hartree - Fock approximation the maximal number of extra electrons < const.) 

Finally we show 

Proposition 1.4. Property (E') implies Property (E). 

Proof. We prove (E) by induction in the number of the atoms k. For k = 2, Property (E) 
follows immediately from Property (E'). We assume that (E) holds for M, replaced by k — 1 
and show it holds for M = k. Indeed, let n\, ...,rik be numbers satisfying the assumptions 
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of condition in (E) for M = k. By relabelling the nuclei, if necessary, we can assume that 
> \ n k\ and n\rik < 0. If n\ > 0, then, by Property (E'), we have that Ei >ni + Ek, nk > 
El,m+n k +Ek. Therefore, I?i jni + ...+i?/ c „i ]nfe _ 1 +£'fc jnfc > Ei :ni+nk +E2 ! n 2 + ---+Ek-i : n k _ l +Ek, 
which together with the induction hypothesis implies Property (E) for M = k. If n\ < 0, 
then, by Property (E'), we have that Ei <ni + Ek iHk > E\p + Ey.^ nk+ni and we arrive to the 
same conclusion. □ 

To simplify exposition we first treat the case without taking statistics into account and 
then modify our results and analysis for the case of an appropriate statistics. Without 
statistics, Condition (D) follows from the positivity improving property of e~^ Ha ,(3 > and 
from Perron- Frobenious theory (see for example [RSIVj ) . Therefore, the statement becomes 

Theorem 1.5 (van der Waals forces without statistics). Consider the hamiltonians H^(y) 
on the entire L 2 (R. 3N ). Then for min{|yj — yA : 1 < i < j < M} large enough (jl.4p holds 
with the constants o~ij defined in (|3.43[) if and only if Property (E) stated above holds. 

Note that when a corresponds to a Young diagram with one row (completely symmetric 
representation) Theorem 11.21 gives Theorem 11.51 

The paper is organized as follows. In Section [2] we discuss preliminaries of quantum many 
body systems. In Section [3] we prove Theorem 11.51 modulo some technical relations which 
are proven in Sections [U [5] and [6j In Section [7] and in Appendix [A] we rework the proof of 
Theorem 11.51 to prove Theorem 11.21 

Notation. We collect here general notation used in this paper. In the text C will denote 
a positive constant which might be different from one equation to the other, but which is 
independent of the nuclear co-ordinates yi, ■ ■ ■ yu- We will use the notation < for inequalities 
that are true up to such a constant. We will write A = B if JB|| < e~ CR , in an appropriate 
norm, where C is some positive constant and 

R = mm{\ yi - y 3 \ : 1 < i,j < M,i + j}. (1.10) 

For any Banach space X, we denote B{X) := {/ : X — > X : f linear and bounded}. 
For an operator A, the symbols o~{A) and o" e ss(^4) stand for the spectrum and the essential 
spectrum, correspondingly. We also use the notation (x) = (1 + \x\ 2 )^ and A = Yl!j=i ^ay> 
V = (V Xl , V^) with A Xj , V Xj the Laplacian and gradient acting on the coordinate 
xj G IR 3 , respectively. Finally, || • || will denote the L 2 — norm of a function or the B(L 2 )— norm 
of an operator, depending on the context, and the symbol 0(5) is understood in this norm. 
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grateful to Volker Bach, Marcel Griesemer, Alessandro Giuliani and Artem Dudko for useful 
discussions. His research was supported in part by NSERC under Grant No. NA7901 and 
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2 Preliminaries about many body systems 
2.1 General properties of Hn(v) 

The general information on the spectrum of the Hamiltonian (|1 .2j) is given in the following 
theorem which is a special case of the HVZ Theorem (see e.g. |HS|, ICFKSj ). 

Theorem 2.1. a ess (HN(y)) = [£, oo), where £ = inf a(HN-i(y))- 

This theorem says that the essential (continuous) spectrum of H^{y) originates from the 
molecule shedding of an electron which moves freely at infinity and therefore whose energy 
spectrum changes continuously. The next result shows that H^{y), as well as each atom, has 
a well-localized ground state (see e.g. [HSl ICSl ICFKS] ) : 

Theorem 2.2. The operator -ffjv(y) has infinite number of eigenvalues, Ej, below its essential 
spectrum, a ess (H]y(y)). Moreover, the corresponding eigenf unctions, ^j, are exponentially 
localized i.e. 

\Vj(x)\<Ce- s \ x \, (2.1) 
for any 5 < */S — Ej . Here x = (x\, . . . , Zjv) an d S := inf a ess {H^{y)). 

The first part of the statement as far as atoms are concerned is known as the Zhislin 
theorem and the second part, as the Combes - Thomas bound. The last theorem shows that 
the atoms and Born-Oppenheimer molecules are stable in the sense that at sufficiently low 
energy they are well localized in the space. However, it says nothing about stability of true 
molecules. The likely source of instability of molecules is not shedding of an electron but 
breaking up into atoms or ions which in total have the same energy as the molecule. So far 
the only molecule proven to be stable is the hydrogen molecule H2 (see |BFGR| ) . 

The uniqueness of the ground state is a delicate issue. Without statistics the ground 
state of Schrodinger operators H is unique (non-degenerate) . This follows from the positivity 
improving property of e~^ H ,/3 > and from Perron-Frobenious theory as mentioned above. 
For spaces with statistics the ground state energy is in general degenerate (for an irreducible 
representation Tg of the permutation Sjv, it is at least the dimension of this representation) 
but its multiplicity is not known. 

Proposition 2.3. If ^ is an eigenfunction of a Hamiltonian Ha (of the form (|2.4p below), 
corresponding to a non degenerate eigenvalue then the one electron density of^ is spherically 
symmetric. 

Proof. For any rotation R in M 3 we consider the transformation Tr defined by 

T R ^( Xl ,...,x k ) = *(i? _1 xi, R~ 1 x k ). 

Since the Coulomb potentials are spherically symmetric, we have that Ha commutes with 
Tr, i.e. HaTr = TrHa- Since * is eigenfunction of Ha the last relation gives that Tr}& 
is also an eigenfunction of Ha corresponding to the same eigenvalue. Since the eigenvalue 
is non degenerate we obtain that Tr^ = c(R)^>, where c(R) is a complex valued function. 
Since Tr is unitary we have that |c(i2)| = 1 for any R and therefore, 

|¥(zi,...,a*)| 2 = \^(R- 1 x 1 ,...,IT 1 x k )\ 2 , 



for any rotation R. Using this and the definition of the one-electron density, we conclude 
that the latter is spherically symmetric. □ 
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Theorems 12.11 and 12.21 still hold for fixed symmetry types, while and an appropriate exten- 
sion of Proposition 12.31 is given in Subsection 17.11 where we begin our analysis hamiltonians 
of fixed symmetry types. 



2.2 Decompositions 

Recall that M and N are the numbers of the nuclei and electrons. Let a = {A\, Am} be a 
partition of {1, 2, N} into disjoint subsets some of which might be empty. With the set Aj 
we associate the j-th nucleus of charge eZj by assigning the electrons of coordinates in Aj to 
be in the same atom/ion as the j-th nucleus. This gives a decomposition of the system. We 
denote the collection of all such decompositions by A and we will call A\,..., Am clusters of 
the decomposition a. The set of all a € A with \Aj \ = Zj for all j = 1, ...,M will be denoted 
by A at . Its elements correspond to decompositions of our system to neutral atoms. 

If a = {Ai, Am} and b = {B\, ...,Bm} are elements of A at , then there exists a permu- 
tation 7T G Sn/S(o) (here / denotes the set theoretic difference), where S(a) is the subgroup 
of <SV consisting of the permutations that keep the clusters A\,...,Am of a invariant, such 
that 

Bj = G A 3 },\Jj € {1,...,M}. (2.2) 

In this case we write b = ira. Various b € A at are related by permutations of the electrons 
and could be labeled as b = a n = ira, tt E Sn/S(cl) with some redundancy. Note that we 
could use also more general K— cluster decompositions and finer decompositions into > M 
subsystems, but we do not need this. For each decomposition a = {Ax, Am} £ A we define 
the Hamiltonian 

M 



H a =J2H Am , (2.3) 



where 



m=l 



- 2 Z ^ e 2 



H Am :=J2(-A Xi -^-) + £ -±— (2.4) 

(so i?A m is the Hamiltonian of the m-th atom or ion and H a is the sum of the Hamiltonians of 
the atoms or ions of the decomposition a) and the inter-cluster interaction I a := H^{y) — H a . 
In other words I a consists of all terms of interaction between the different atoms/ions in the 
decomposition a. We have that 

H N (y) = H a + I a . (2.5) 

Let E a and ^ a be the ground state energy and ground state of H a , respectively. Then 
H a ^ a = E a V a . We have that E a = Y^ =l E Am and 

M 

^ a (x x ,...,x N ) = (pA m (xA m ), (2.6) 

m=l 

where E& m and <j>A m are the ground state energy and the corresponding ground state of the 
atom/ion A m and XA m = (xi : i € A m ). (For an atom we denote -EU m ='■ E m .) Furthermore, 



4>A m {x Am ) = (t>m{xA m ~ Vm), 



(2.7) 
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with (j) m the ground state of the m-th atom with the nucleus fixed at the origin and XA k —y m = 
(%i — Urn '■ i £ A m ). Throughout the text we will always assume that ||^ m || = 1 for all 
m = 1,...,M. Note that by the definition of a's, min ag ^ E a = E{oo), where E{oo) was 
defined in (jl.3p . (If Property (E) holds, as we expect it always does, than E a = E(oo), 
Va 6 A at .) Standard estimates (see |HSj ) give that there exists > such that 

\\e d{xA ^d a (p m \\ < l,Va with < \a\ < 2, (2.8) 

where a is a multiindex with each index corresponding to differentiation in some body vari- 
able. For a = {Ai, ...,A M } G A we define y a = (y a>1 , ...,y a ,N), where y a>i = y m for i G A m . 
In other words y a ^ is the coordinate of the nucleus that Xi is assigned to. From (|2.6p . (|2.7p 
and p.8p we obtain that 

\\e e{x - ya) d a ^ a \\ < l,Va with < |a| < 2, (2.9) 

where x = (x±, xjy). 



3 The general approach (without statistics) 

In this section we lay out the general set-up for our approach without taking the statistics 
into account. Our goal is to prove Theorem 11.51 under certain technical assumptions, which 
are then verified in later sections. In what follows we omit the argument y and write E and 
H for E(y) and Hiv(y), respectively. 

3.1 Feshbach map 

Let P be an orthogonal projection and P 1 - = 1 — P. Introduce the notation H = P^HP^. 
We will use the Feshbach-Schur method (see [BFSt IGS] ). which, as applied to the quantum 
Hamiltonian H, states that if 

(a) Ran(P) C D(H) (domain of H) and therefore ||-ff-P|| < oo; 

(b) The operator (H 1 - — A) is invertible; 
then the Feshbach-Schur map 

F P (\) = (PHP - [/(A))| RanP , (3.1) 

where 

U(X) := PHP ± (H ± - \y 1 P ± HP, (3.2) 

is well defined and 

A eigenvalue of H <^=> A eigenvalue of Fp(\). (3.3) 

Moreover, the eigenfunctions of H and -Fp(A) corresponding to the eigenvalue A are connected 
as 

Hip = \ip & F P (X)(j) = \<f>, (3.4) 
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where (j>, ip are related by the following equations 

(j) = Pip, V = Q(X)<f>. (3.5) 

Here the family of operators Q(X) is defined as 

Q(A) = P - P ± (H ± - X)- 1 P ± HP. (3.6) 

(Remember that we do not display the y— dependence of various objects.) 

In what follows we will explain how we are going to use the Feshbach map. The first step 
is to choose the orthogonal projection P. 



3.2 Orthogonal projection P 

We define the projection P, to be used in the Feshbach-Schur method described above, as 
the orthogonal projection on spanj^ : a G A at }. According to (|3.ip . we have to compute 
PHP and U(X). To this end we use that the functions ^f a are almost pairwise orthogonal 
for large R (where R is defined in (jl.lOp ). i.e that 

(Va,Vb} = Sa,b, (3-7) 

following from (|2.9p . This implies that the matrix {g a b) '■= ((^a> ^b)) = 1- Hence, we can 
write 

P= iW^fcl, (3.8) 

a,b£A at 

where (g ab ) is the inverse of the matrix (g a b)- Since (g a b) = 1 5 we have that this inverse 
satisfies 

(g ab ) = l. (3.9) 

As a consequence, we have 

P= l*»X*«l- ( 3 - 10 ) 

a£A at 

To show that the Feshbach map exists, we note that Ran(P) C T)om(H) since the range of 
P is spanned by {^ a ) a eA at an d by (|2.9p each ^ a is in H 2 (M. 3N ). Hence the condition (a) 
holds. We will prove in Section [J] that, under Property (E), there is a 7 > 0, independent of 
y s.t. for R large enough, the stability bound 

H 1 - > E(oo) + 27, (3.11) 

where E(oo) was defined in (|1.3p . holds. We assume this for now. Using the last relation we 
obtain that H 1 - — A is invertible for all A < E(po) +7. Thus the condition (b) is also satisfied 
and, for all A < E(oo) + 7, the Feshbach Schur map -Fp(A) is well defined. 

Our goal now is to estimate the two terms on the r.h.s of (I3.ip . Before we proceed we 
record the important fact that P commutes with the permutations 

PT n = T n P, Vvr G S N , (3.12) 

where, recall, are unitary operators given by (11.51) . 
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3.3 Estimate of PHP 

We show that 

PHP = E(oo)P. (3.13) 
Indeed, using f|3.8[) and (|3.9p we obtain that 



PHP= \^a}9 ab (^b,H^ c )g cd ^ d \ 

a,b,c,d£A at 

= ]T \* a ){*a,H* c ){* c \. (3-14) 
a,cg.A at 

Next, the relations H^ c = (H c + 7 C )* C = (-E(oo) + 4)* C! ([377]) and (f3~T0|) give 

PHP = E(^)P+ Yl l*»X*a,W<*c|- (3-15) 

a,cGA at 

Lemma 3.1. For all a, b € A at , we have 

(* a ,I b * b ) = 0. (3.16) 

Proof. We first prove (|3.16j) for a ^ b. By Cauchy Schwartz and the inequality \\fg\\ < 
||/||loo ||<y|| (recall that for functions |.| is always the L 2 norm) we obtain that 

Now, using (|2.9p to bound the first factor on the r.h.s., using the estimate e~ 9 ^ x ~ Va ^ +( ' x ~ yb ^ < 
e -d(y a -Ub) ; which follows from the triangle inequality for the bracket (•), to bound the second 
factor and finally using that If, is bounded relative to the Laplacian to handle the third one, 
we obtain that 

K*«, Wl < e~ cR ||(l - A)e^-^\P 6 ||. 

The last relation together with relation (|2.9p and the boundedness of the derivatives of (x—yb) 
implies relation (13.16P for a ^ b. 

We will now prove that = 0. Let a = {A%, Am}- We have that 

x,M 

where 



e 2 e 2 e 2 e 2 



4" = -rm ~ — r + r~ i + — r ■ ( 3 - 18 ) 

Wj xi\ \yi x m \ \tji y,j\ \x\ x m \ 

We pass to the variables 

zi r = xi-y r , VI 6 A r ,Vr = 1, ...,M, (3.19) 
and write za t = (zi r : I € A r ) and dzA r = YlieA r dzi r . Using (|2.7p . we obtain that 

(Vajfr**) = / dz^dZA^f |^(^)| 2 |^(^)| 2 , (3.20) 



Long Range Behaviour of van der Waals force, April 22, 2012 



14 



where 

=2 2 Jl Jl 



j lm = e f |_ _f i f ; (3-21) 

\Vij + Zli\ \Vij ~ Z m j\ \Uij\ \Vij + Zli Z m j | 



Since by Proposition 12.31 the one electron densities of the functions <f>i , <f>j are spherically 
symmetric, we have by Newton's Theorem (see for example [LL] ) that 

H^)| 2 1 dZAi =J-rf \<P(ZAj\ 2 dZ Al . 

\Ulj -T z h\ \yij\ J \g u \<\y i: j\ 

with yij = yi — Uj. Using that the ground state is normalized and is exponentially decaying, 
we have 



/ l^(^)| 2 | l ——Az Ai = ^- (3.22) 

J \yn + zii\ \yij\ 



\Vij + zu\ \yij 

As a consequence, 

1 1 



|^(^)| 2 |^(^.)| 2 ^ — dzAi dz A . = ^,VZ G Ai. 



In a similar way we can obtain that 

1 

I Z m j Vij | | yij I 

Furthermore, using (|3.22[) twice, we can show that 



l^^.)! 2 !^-^.)! 2 - -,dz A .dz Aj = - — r ,Vm € Aj. (3.24) 



\(t>i{z Al )(t>j{z Aj )\ 2 1 _ -dz Ai dz A . = j—. (3.25) 

I y^j ' Zli Z m j | | yij I 

Equations (|330|) . ^21} . (I^23|> . (jEMj) and (I3T251) imply (tf 0) ig^a) = 0, for all i,j G 
{1, ...,M} with z ^ j and / G Aj,m G Aj, which together with (|3.17[) implies that (^ a , I a ^ a ), 
and therefore (I3.16|) is proven. □ 

The relations (|3351) - (I3T61) imply (jBTTHjl . 
3.4 Rough bounds on E 

Before proceeding to analysis of U (A) we prove rough bounds on E. 

Lemma 3.2. The ground state energy E of H satisfies the following inequalities: 

M 

x ±^\Vi-Vi\* 

Note that the left hand side gives a bound from below of the van der Waals interaction energy 
which is sharp in the order of magnitude. 
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Proof. The right hand side of (|3.26[) follows immediately, since using the relations (|2.5p . 
H a ^ a = E(oo)^ a , and (13.161) we obtain that 

E < (¥ , HV a ) = (¥ a , P a * a ) + / a * ) 

= P(oo) + (* a , J ¥ a ) < E(oo) + e" * 

We will now show the lower bound in (|3.26p . 

Let A < E(oo) + 7, where 7 > is the same as in (pTTT]) . Then by and (pTTTj) . [7(A) 
is well-defined and is positive definite, U(X) > 0. Therefore, the finite rank operator Pp(A) 
defined in (|3.ip is well defined. Let v{\,y) be the lowest eigenvalue of Pp(A) (recall that H 
and therefore Pp(A) depends on y). We consider the fixed point problem 

A = v(\,y) 

on the interval (— oo,P(oo) +7]. We will show below that for R large enough the fixed point 
problem has a unique solution Ao(y). 

Using j£g)> CT, = CHa + 4)*a = (P(oo) + and that P^oo)^ = 0, we 

obtain 

P X PP = P L \Ia*a)(*a\ and PPP 1 = | * ) (J * a |P X . (3.27) 

(Note that due to (|2.9p the operators on the r.h.s and l.h.s. are bounded.) The relations 
(pTZD , (|3TT]1 and (j33Tjl give that 



a,b£A 

where 



C/(A) = ^ |* a )C/ o6 (A)(^ 6 |, (3.28) 

a,beA at 

U ab (\) := (* a , I a P L {H L - A)- 1 P ± / 6 * 6 ). (3.29) 



We show now that 

l,M 



Recall that for a = {A\, Am}, the intercluster interaction I a can be written as (|3,17p 
- i.e. J a = E^^, ^ = Eie^meA^lT' with 4" S iven in The 

function /'j™, written in the variables (|3.19p . is 

7T™ = 1 1 . 

I Zli Uij I I Uij ~ z mj | | Uij \ \ z li z mj Uij \ 

Using that 

| Z± ,«|-' = fe|- 1 d ± ^f + ^)-, (3.31) 
and expanding in zn and z m j, we obtain that 



p 2 

jlm ° 

ij 



\y%j\ 



(zu z m j) • yij zn • yij z m j • yij \zu z m j\ 3((zn z m j*) • yij) 



21 



\Vij\ \Vij\ \Vij\ 2 \Vij\' 



Long Range Behaviour of van der Waals force, April 22, 2012 



16 



,2 



\Vij\ 



\zii\ 3(zn ■ Uij) \Zmj\ 3(z m j ■ Uij) 



+ 



2\Vij\ 2 2\Vij 



2 



+ 0{ 



\z\ 3 



\Vij\ 



where \z\ := max; im (z;j, z m j), provided \z\ < k\Uij\, which in this case gives that 

4 m = ltW?^ ^ + °(t^t4)' for N^b*;l> ( 3 - 32 ) 

IVijl IVijl 

where f-J 1 := zu ■ z m j — 3(zu ■ yTj)(z m j ■ y~ij). On the other hand, using (|2.9p and using that 
I a is bounded relative to the Laplacian, we obtain that 

Il4x k |>|| 2/i .|^a||<e- c ^ 1 , (3.33) 

where xci is a smoothed out characteristic function of the set 0. This estimate together with 
(13321) gives (13301) . The estimate (pTSPj) . together with (03} , (f3T28l) . (13391) . (((^ - A)" 1 )] < 1 
and the arithmetic mean-geometric mean inequality, shows that under (|3,llj) . 

M 

\d?U(\)\ < ^ — m > 0, (3.34) 
, ■ . ■ \yi yj\ 

uniformly in A < E(po) + 7. 

By iH), iHD and flffl), we have that |i/(A, y)— E(oo)| < Eiii<i l^jpr and |0 A f (A,j/)| < 
2i=i<j |y. .|6 uniformly in A < E{oo) + 7. Hence, for i? sufficiently large, z/ is a con- 
traction and maps (— 00, i£(oo) + 7) to \E(pa) — j,E(oo) + 7]. Since v maps the interval 
[E(oo) — 7, E(oo) + 7] to itself and is a contraction on this interval, the fixed point problem 
has a unique solution Xo(y) in [E(oo) — 7, E(oo) + 7]. 

From IpTL]) . (|3"7i"3j) and (f334|) we obtain that 

^ 1 

F P (Ao(y))-£(oo)> 



and since Ao(y) is eigenvalue of Fp(Ao(y)) we obtain that 



M 

;l2/i-%| 6 

By isospectrality of the Feshbach map, (13. 4j) . we obtain that Ao(y) is an eigenvalue of H. 
We will now show that 

E = A (y), (3.36) 

i.e that Ao(y) is the ground state energy energy of H. Indeed, since Ao(y) is eigenvalue of H 
and E is the lowest eigenvalue of H we have that E < Xo(y). So it is enough to show that E > 
Ao(y). First we show that for fixed y, Fp(\) is a decreasing function of A on (—00, E(oo) + 7]. 
Indeed, for any Ai,A2 € (— 00, .E(oo) +7] with Ai < A2 the operator (H 1 - — A/)" 1 is defined 
bounded and differentiable on [Ax, A2] and therefore by the Fundamental Theorem of Calculus 



A' 

2 



(H x - A2)- 1 = (H x - AO" 1 + I \h x - A) 

JAi 



dx. 
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Since (H 1 - — A)~ 2 > 0, we obtain that (H 1 - — A/)" 1 is an increasing function of A on 
(— oo,i?(oo) +7]- Therefore, by (|3,ip and (|3.2I) Fp(X) is a decreasing function of A. Since 
E < Ao(y) and -Fp(A) is decreasing in A we obtain that 

F P (X (y)) < F P {E). (3.37) 

Since we have Fp(Ao(y)) > Ao(y) (because Ao(y) is the lowest eigenvalue of Fp(Xo(y))), the 
relation (|3.37p gives that Xo(y) < Fp{E). The last relation together with the fact that E is 
eigenvalue of Fp(E) gives E > Xo(y). The relations E < Ao(y) and E 1 > Ao(y) give (|3.36j) . 
The relations f)3.35j) and f|3.36j) give the left inequality of (|3.26j) . This completes the proof of 
(PL26D - □ 

The estimate (|3. 11)) together with (13.261) gives that for R sufficiently large 

H x - E > c> 0, || {H L - E)- 1 1| < 1. (3.38) 

Consequently, due to its definition (|3.29j) . it follows that the matrix (U ab (E)) is well-defined, 
non-negative and uniformly bounded. 



3.5 Estimate of U(E) 

We recall that H^ m denotes the Hamiltonian of the m-th atom and E m denotes its ground 
state energy. We let Py stand for the projection onto the ground state of Ha. + Ha } and 

p ij '■= 1 - Piy Let l^tjl = \Vi ~ VjV 

Recall that x = (xi, ...,xn) is the collection of the electron coordinates. Let <p(x) = ^r{x) 

be a C 2 monotonically increasing function of (x), where, recall, (x) = (1 + |x| 2 )2, with 

uniformly bounded derivatives of first and second order (in both x, R) and satisfying, 

, . / (x), if (x) < R - 1 , A 

<Pto = { R, if(x)>R + l. (3 - 39) 

Recall that y a = (y a ,i, y a ,N) where y a ,i = y m for i £ A m and a decomposition a = 
{Ai, Am} G A at . In other words y a ^ denotes the coordinate of the nucleus that the 
electron with coordinate is assigned to by the decomposition a. For any 6 > 0, we let 

We show in Section [5] that E is in the resolvent set of H^ and 

\\{Hi-E)- l \\ <1. (3.41) 

Our goal is to show 

Lemma 3.3. Assume (|3.4ip holds. Then there are positive constants, independent of y 
and depending only on the nature of the atoms i and j, s.t. 

l,N 4 

U(E) = f(y)P, where f(y) = £ ^(1 + O(-)). (3.42) 

2<J 
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Moreover, are given by relations 

o-ij = (fijfcfa, (P t j(H Ai + Ha.)P± -Ei- Ej)- 1 firfrfj) (3.43) 
where fa := Eke^leA, /§'• 

Proof. We use (pPdjl to estimate £/ ab (£) in (|3T29]l . We claim that 

U ah {E) = U aa (E)5 ab (E) and independent of a Va,6e„4 a *. (3.44) 
Indeed, using (|3.8p . (|3.9p and (|3.16p . we obtain that 

P-L/a^o = ^a*a, ^6*6 = V^6- (3.45) 

This, together with (13. lip . (I3.29P and the relation 



e 



-6<p(x-y a )( H ± _ E yl e S V (x-Va) = ( H ± _ E y^ ( 346 ) 



where 5 > 0, implies that we have that 

U ab {E) = (J tt tt 0> (H L - E)- 1 ^,) (3.47) 
= (eM'-^IaVa, (Hi- - E^e-^-^I^b). (3.48) 

Due to (|2.9p . and the Laplacian boundedness of I a we have for all 5 < 6 that 

|| e M*-*0j tt tt o || < 1 Va. (3.49) 

In addition, by the construction of (p (see (I3.39P ) we have that ip(x—yk) > min{(x— yk), R— 1}. 
Using this and the triangle inequality (x — y a ) + (x — y&) > (y a — y^) > R (since a ^ 6), we 
obtain the inequality ip(x — y a ) + tp(x — yb) > R — 1, which in turn implies 

e -6<p(x-y a ) < e -8R e Sip(x-y b ) ^ (3.50) 

Using the relations (|3ldl . (EOS]) . (I3^9|) and (1330]) we arrive at (13131) Va ^ 6. 

Now we show the a = b part of (|3.44p . i.e. the independence of a. First, recall that for 
any a,b £ A at there exists a permutation 7r with b = ira (this equality was defined precisely 
in (HOI) ). Moreover, = vE^. The last two facts imply that 

T v I a * a = J 6 * 6 . 

Using (|3.29p . and the last relation we obtain that 

U bb (E) = (T„I a V a ,P- L (H J - -EY l P^I a ^ a ). 

Since H and P commute with T n (see f)l . 7f) and (|3.12p ) and since T n is unitary, the last 
relation implies that U bb (E) = U aa (E), i.e. the a = b part of (EOip . 

Thus it suffices to estimate U aa (E) for a fixed a. Recall the formula (13.17P (I a = 
^2\Li <m lAiAmi where lA t A m are the interactions between the atoms Ai and Aj). The es- 
timates (|3.32p and (I3.33P imply that the terms I Ai Aji entering (|3.17p . satisfy the estimate 

I Ai A^a = j^Mz^a + O^), (3.51) 
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</?!;$ 4 if « 7^/3, (3-55) 



where the remainder here is in the sense of the I? norm, y^j = t^t, z := (zu, z m j, I € A{, m € 
Aj) with the variables zn,z m j, defined in (|3.19p . and 

fij(z, Vij) = ^2 {zki ■ zij - 3(z ki ■ yij)(zij ■ yij)). (3.52) 
We insert (|3"T7|) with (|33Tj) into (13^47]) with b = a and use (13381) to obtain 

^)-E^^(E^). ("3) 

where a, /3 run in pairs of nuclei in the decomposition a and 

<p = (f a *a, (H X - EY l f^ a ). (3.54) 

We prove in Section that 

\rr a n\ < 

R 

and that 

a - a = aa + (^), (3.56) 

where a a , for a = (ij), are given in (|3,43p . 

The relations (E2SD, ([531D , (13331) . (|335j) and ([336]), together with (13TTU]) . imply <pP2|> 
modulo the fact that o"jj are positive and independent of y which we prove below. 

First we note that the only dependence of au on£ := , v ' V] appears on so we will 

write fjj j and . For any rotation R in R 3 we define Tr acting on the space L 2 (Mp^ Ai ^ Aj ^) 
(recall that \Ak\ is the number of electrons of the &;-th atom) of functions of the variables 
Zki = Xk — yi, k = 1, . . . , \Ai\ + \Aj\, I = 1,2, by rotating them. Since Tr commutes with 
i//^ + Ha 3 and since is a ground state of Hj^ + Ha^ , we have that Tjifacfrj is also a 
ground state of il^ + Haj • Therefore, since the ground state is non degenerate and Tr is 
unitary we obtain that TRfafy = c{R)<f>i(j)j where \c(R)\ = 1. This implies 

Tr-P^- T R = . 

Using the last relation and the fact that Hj^ + commutes with Tr, together with (|3.43p 
and the fact that Tr is unitary, we obtain that 

of = Vr Af^CO^.il'jilh, I H A .)P+ I-, l-p '//, Ajfj ^o,)'. (3.57) 
On the other hand using (|3.52p and the fact that R is unitary we obtain that 

Vr'./;' J'u "'• (3-58) 

But using that T^fafa = ciR' 1 )^, with \c{R^)\ = 1, together with fl3J)8j) ; we obtain 
that 

Tariff 1 Mi) = <R 'If!; ;?, '"-Or 
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The last relation together with (|3.57j) and \c(R 1 )\ = 1 implies that 

4' = 4 '*'OiO,, (/',;• (//.i. + HaM'j -Ek- E^uf^Mi)) = of\ 

which implies that Oij is independent of g/j,-. 

Thus, it remains to prove that o~ij > 0. By (|3.43p and the fact that PA(E[ Ai + Ha^P^ — 
Ei — Ej is a positive operator, we obtain that o~ij > 0. Assume now that o~ij = 0. Then 
since (P-j^HAi + Ea^)P^ — Ei — Ej)" 1 > 0, we obtain that fij<j>i4>j is a minimizer of the 
quadratic form ('J/, {P^-{R~a 1 + Ha^P^ — Ei — Ej) -1 ^) and therefore the Gateaux derivative 
vanishes at fijfafa or in other words {P^{Ha. + Ha^P^ -Ei- Ej)' 1 fijfafy = 0. Applying 
(P^-(Ha 1 + HAj)Pjj — Ei — Ej) on both sides of the last relation we obtain that fij<fii(f)j = 0, 
which contradicts the fact that fijfafy ^ 0. □ 



3.6 Conclusion of the argument 

From (j3TT]h (pnSll and (l3^2|) we have that 

F P (E) = E(oo) - f(y). (3.59) 

Since the ground state energy E is an eigenvalue of H, we conclude from (|3.3p . (|3.59p and 
the definition W(y) = E — E(oo) that (11.4j) holds. Therefore, we have proven the following 
result: 

Theorem 3.4. Assume Property (E). If the relations ([37TT]) . (j3Hj) . (j335j) and (j336l) /ioZd, 
i/ien so does the van der Waals law without statistics (Theorem \1.5\) . 

3.7 Proof of the necessity of Property (E) 

In this section we will show that if Property (E) fails to hold then so does the van der Waals - 
London law. We denote by .4 mm the set of all a € A for which inf cr(H a ) = E m i n . If Property 
(E) holds if and only if A mm = A at . We now use the Feshbach map but with the orthogonal 
projection P onto 

span{^ a : a € .4 min }. 
Then proceeding as in the proof of (|3.10p we can show that 

P= Yl l*a)<*o|- (3-60) 
aG.A min 

Furthermore, repeating the arguments of the proof of (|3.1ip (with A at replacing >l mm ) we 
can prove that there exists 7 > such that H ^ > E m i n + 27, which implies that the Feshbach 
Schur method can be used. 

We will now estimate PHP. Using (13.60P and proceeding as in the proof of (I3.15p . we 
can show that 

PHP = E min P+ l*a)(*o,W(*c|. 
a,cG-4 min 
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Due to exponential decay of the ground states and the result (|3.16p , we have that 

(* a , I c * c ) = when either a / c or a = c E A at . 

Since the Property (E) is not satisfied it means that there exists a £ A mm with a ^ A at . For 
such a decomposition a = . . . ,Im), we define the charges q a i := (Zi — |ij|)e. Next, we 
apply Newton's theorem (13.22h to show that 

Clearly, ^ j^y x Ti ■ The last three relations imply that 

PHP = E min P+ V r^W^I- (3.61) 

Now, we estimate U(X). Proceeding as in the proof of ([3.28P we can obtain that 

17(A) = \y a )(Ia*a,P ± (H ± -\r 1 P ± I b * b )(y b \. 

a,&€.4 min 

In addition proceeding as in proving (13.511) . we can show that 

\\I*Va\\<^,VaeA, (3.62) 
The last two relations together with i7~ L > E m i n + 27, for 7 > 0, give 

U(X) = 0{jp), VA < E mhl + 7. (3.63) 
The relations ([33]) . ([33T]) and (1333]) give that 



F P (X) = E min P+ V V^^|^a)(^a| + 0(4)' VA<75 min + 



7, 



a,fce.4 min /-4 at i^j 



and therefore, by Yli^j \y % ~ a y- \ x \ anc ^ the van der Waals law fails. As a consequence 

Property (E) is necessary for van der Waals law to hold. 



4 Proof of stability bound (13.111) assuming Property (E) 

Recall that the system we consider consists of M atoms with Z\, Zm electrons, respectively, 
and that E m and E m ^ n denote the ground state energy of the atom m (corresponding to the 
m-th nucleus) and the ground state energy (the infimum of the spectrum) of the ion with a 
nucleus of charge eZ m and Z m — n electrons, respectively. Let Ea be the first excited state 
energy of the atom j and E(oo) as defined in (jl .3[) . We define 



71 = min{£'. - Ej : j = 1, M}, 



(4.1) 
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and 

72 = min{(£ 1>m + ... + E M ,n M ) - E{oo)}, (4.2) 

where the minimum is taken among all n\, ...,um £ Z, satisfying the constraints of Property 
(E). Due to (|1.3p and Property (E) we have that 72 > 0. Finally, we let 

70 := min{7i,7 2 } > 0. (4.3) 

The relation (|3.1ip will follow immediately from the following proposition: 

Proposition 4.1. H 1 - > E(oo) + 70 — o~(R), where S(R) — > as R — > 00. In particular, 
(|3,lip holds for R large enough. 

Proof. Recall that A is the set of all decompositions, which correspond to all partitions 
a = {Ai, Am} of {1,...,N} with \Aj\ > for all j. Thus, A includes decompositions of 
the system with clusters that are possibly ions. We will split the configuration space R* N 
into domains. For each a = {A±, Am } € A, we define 

fi£ = {( Xl ,...,x N ) : \xi- yj \ > (3R,Vj = 1,... I M,V» Aj}. (4.4) 



Following Sig we will now construct a partition of unity ( J a )aeA having the properties: 

< J a < l,supp J a C ttl,\\d a J a \\ L °o < R-W, (4.5) 

J a commute with all elements of S(a), (4-6) 

where, recall, S(a) is the group of permutations that keeps the clusters of a invariant, for all 
a € A, and 

E J « = L ^ 

aeA 



We consider the functions F a = x 3 * <fi where 4> : R — > ]R is a function supported in 

Bj_(0) (ball of radius centered at the origin) symmetric with respect to particle coordi- 
nates, with (p > and J R3N <p = 1 and xa denotes the characteristic function of the set A. 
Then F a € C°° and F a > 0. Furthermore using the triangle inequality and the fact that 4> is 

supported in Bj_(0) we obtain that supp(F a ) C f2a and that F a \ 2=1. The last relation 

10 n| 
2 

together with the fact that Ua^A^a = M 3Ar gives that ^ a ^F a > 1. Therefore, if we define 

j F « 

'J a 



then J a is a partition of unity satisfying all the desired properties. 
Now we use the IMS localization formula (see for example |CFKSj ) 



H = Y,{JaHJ a -\VJ a \ 2 ). (4.8) 

a£A 
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This together with (|4.5p gives that H > X^agyl JaHJ a — %■ The last relation together with 
([23]) and I a > - Ysf=i Sili.igAj - ~% on sn PP J a implies that 

H>J2 J aH a J a -^- (4.9) 



The last relation together with the fact that P is self-adjoint and P < 1 implies that 

C 
R 



P ± HP ± > £ P^JaHaJaP 1 - ~. (4.10) 



We will now estimate each of the terms on the right hand side of (|4.10p . 
Lemma 4.2. With 70 defined in (|4.3p . we have that, Va E A, 



P^JaHaJaP^ > {E(oo) + 70 )P ± J^P 1 - ^. (4.11) 

-ft 

Proof. We will consider the following two cases 

Case (1): a € *4/„4 a * (i.e. some clusters of the decomposition a include ions). 

In this case we define rij := Zj — \Aj\. Since by definition Ej >n . is the infimum of the 
spectrum of the hamiltonian -Ha,, defined in (|2.4p . we have that > Ej jnj . The last 

relation together with (|4.2p . (|4.3p and (|2.3p gives that 

#a > + ... + E M ,n M > (E(OO) + 70 ) 

and as a consequence 

P L J a H a J a P L > (E(oo) + lQ )P L J 2 a P L . (4.12) 
Case (2): a € A at (i.e. a is a decomposition of the system to M (neutral) atoms). Now, let 

Pa = |*o)(*a|, (4-13) 

where recall that ty a is the ground state of H a . Using (P 1 - — P^)J a H a = (P a — P)J a H a , 
(|3.8p . (|4.13p and the self-adjointness of J a and H a , we obtain that 



{P ± -P^)J a H a = \V a )(H a J a y a \- \^d)g db (H a J a y b \. (4.14) 

d,b£A at 

We will now prove that 

\\H a J a * 6 || = 0, Va, 6 € „4 a *, a ^ 6. (4.15) 
Indeed, by the Laplacian boundedness of the Coulomb potential we have that 

\\H f Mb\\ < \\J a *b\\ + ||A(J a * 6 )||. (4.16) 

1 

Since J a is supported on Q% , by (|4.4p and (|4.5p we obtain that 

|5°J a |e-^-^ < |9 Q J a |e-¥ < e -¥, (4.17) 
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where 6 is the same as in (|2.9[) . The relations flUTTj) and (|2.9|) imply that 

\\d a J a d^ b \\ = o. 

which together with (|4.16p implies (|4.15p . The latter inequality and (|4.14p yield 

W^P 1 - l£)J a H a \\ =0. (4.18) 
Similarly we show that \\H a J a (P — P^)\\ = 0, which gives 

P^JaHaJaP 1 - = P^J a H a J a P^. (4.19) 

Next, using the relation P^J a = JaPa + {P<ti Ja] = JaPa ~ [Pa, Ja], we obtain that 

Pa JaHaJaPa — JaPa ^-oPa Jo, [*^o> Pa]HaJaPa JaPa J^a[Pa, Ja]- 

Since P a = \^ a ){^ a \ and H a is self-adjoint we have that 

[J a ,P a )H a = \J a V a )(H a V a \ - \V a )(H a J a V a \. 

Using now that H a ^ a = E(oo)^ a and H a J a ^ a = E(oo)J a ^ a + [H a , J a ]^ a and the relation 
(14. 51) . we obtain that 

[J a ,P a ]H a = E(00){\(J a - l)V a }(*a\ " \*a}((Ja ~ l)*a|) + O(^). 

(Our definition of J a (see the paragraph after (|4.7|) ) implies that |<9 Q \l/ a | < e~ cR on supp VJ a 
so that the error above is in fact 0{e~ cR )). By (j2.9|) . we obtain that 

||(1 - J a )* a || = ||(1 - J a )e- e ^-y^\\ L oo\\e 9 ^-y^ a \\ < ||(1 - J a )e^ x ~y^ \\ L «,. 



Using the relations < 1 - J a < 1 - Jl and (@2D, we find that (1 - J a ) < Y,baA,b^a J b- Now > 

R 

5 



since by (|4.5j) . (x — y a ) > -P on the support of Jb for all b ^ a, we obtain that 



0< Jb^ e{x ~ Va) < E Jbe-°r< e - 9 T, (4.20) 

where in the last step we used (|4.7|) . The last three relations give 

||(1- J„)*„||=0. (4.21) 

This shows 

\\[Ja,Pa]H a \\<^. (4.22) 

The relation ||iT a [P a , J a ]\\ < can be shown similarly, which implies 

P^JaHaJaPt > J a P^H a P^J a - |. (4.23) 

Now, if a = {Ai, Am} then by (|2.6p we can write P a = YifLi^Aj where Pa. is the 
orthogonal projection onto the function defined in (12.7p . Since for all j = 1, M we have 
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that P A . commutes with H a and P^.P Aj = P Aj P^ = and P± = E^Pj. 11^+1^)- 
we obtain that 

M M M 

ptH a pt =E( p i( n p **) h j%( n p ^))- 

j=l i=j+l i=j+l 

The last relation together with the relation (|2.3|) for the decomposition a gives that 

M M M M 

ptH a p^ =£E(( n ^j^t^pi-^ n ^))- 
i=i i=i »=j+i i=j+i 

Since all the projections on the right hand side commute and H Al > E\ (recall that E\ is the 
ground state energy of the atom corresponding to the nucleus I) and P^ Ha^a > {Ei+lijP^ 
(by the definition of 71 in (|4.ip and the fact that P^ removes the ground state energy from 
the spectrum of Ha { ), we obtain that 

M M 

PtHaPt > E( II Pa 1 )Pa j {E{^)+ 1i ) = P^(E(oc)+ 7l ). 

j = l i=j + l 

The last relation together with (|4.3|) completes the proof of 

J a PtH a P^J a > (E(oo) + l Q )JaPaJa. (4.24) 

Using that {Pa) 2 = Pa an d that [P^, J a ] = —[P a , Ja) we obtain that 

JaP a Ja = Ja(P a ) Ja = P a JaPa + [Pai Ja]P a Ja ~ P a Ja\Pai Ja] 

The last relation together \\[P a , J a ]\\ < j=r, which can be proven similarly to (I4.22j) . give that 

J a Pa L Ja<PtjlPa L + ^ (4-25) 

Proceeding as in the proof of (|4.18j) we can obtain that 

||(P - P a ) J a \\, || J a (P - P a )\\ = 0. (4.26) 
The relations P~25]) and (f¥726|) imply 

J a P a L Ja<P ± JiP ± + ^. (4.27) 
Finally by ([Q]> . (j4~T|) and ([43]) we have that 

P(oo) + 70 < #1 + P2 + ••• + E M - 

Using the last relation and the fact that the ground state and first excited state of every 
atom is always negative we obtain 

£(oo)+7o<0. (4.28) 

Equations (|4TT3j) . (ET23D . (gjZjJ) , (|4 ^27D and (ET281) imply (|4~TTjl . This completes the proof of 
Lemma 1131 □ 



Equations P~7]) . (14~TU1) and (fiTTl) give that H 1 > (E(oo) + 70^ - %, which together 
with P 1 - < 1 and (|4.28p gives that H 1 - > E(oo) + 7o — This completes the proof of 
Proposition 14.11 and therefore of (|3.1ip . □ 
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5 Proof of ( ED) 



The proof consists of two lemmas. Recall that E = E{y) is the ground state energy of H and 
A = J2j=i A x r We will often use the bound 

||(l-A)( J ff # + C)- 1 ||<l, (5.1) 

where H# is one of the Schrodinger operators which appear in the proof and C is such that 
H# + C > 1. This bound follows from the fact that the Coulomb potential is bounded relative 
to Laplacian with the relative bound zero. For any operator K and 5 > 0, we let 

Kg := e -6 V (x-y a ) Ke 6 V (x-y a )^ R ± ._ ( 5 . 2 ) 

Lemma 5.1. We have that 

WiHs-HXl-Ay'W^S. (5.3) 

Proof. Observe that 

H s -H=-A s + A. (5.4) 
since A 5 = e^Ae - ^, where <p was defined in (13.39[) . by the Leibnitz rule we obtain that 

-A s + A = 6[(A<p) + (Vp) • V - S\V(p\ 2 ]. 

Since ip by definition has L°° bounded derivatives this implies that [(Aip) + (Vy>) • V — 5\ V<^| 2 ] 
is —A bounded. Hence, we obtain (15.31). □ 



Lemma 5.2. The following relations hold for small 5: 

\\Ps-P\\<5 (5-5) 

and 

\\H S (P S -P)\\<6. (5.6) 
Proof. Since, by lETBjl . P = J2b,ceA at I ^ r &>£' 6c <^ r c | , we have that 

p s -p= E (m}9 bc (Vc\)6-m9 bc (y c \)- (5.7) 

b,c€A at 

We show that 

\\(\y b )g bc (y c \) s -\V b )g bc (y c \\\<6. (5.8) 
We start first with the case b ^ c. Let 

g {5) = e - s ^ x ~ y ^\^ b )g bc (^ c \e s ^ x - y ^ - \^ b )g bc {^ c \. 

Then g(0) = and 

g'{5) = -<p{x - y a )e- s ^-y*)\i£ b )g bc (* c \eW*-^ 
+e-W x -y^\* b )g bc {*Mx - yJeM*-*^ 



Long Range Behaviour of van der Waals force, April 22, 2012 



27 



Using now that < tp(x — y a ) < R and \g bc \ < e ( CR > (because of (|3.9[) and the fact that 
b ^ c), we obtain that 

\\g'{8)\\<2Re^ R \\^ b \\\\^cl 

As a consequence for 5 small enough we obtain that ||<7'(<5)|| < 1 and since g(0) = we obtain 
that < 5. This shows (|5.8p for b ^ c. We now consider the case b = c. 

Since g bb = 1 it is enough to show that 

\\(\* b )(* b \) s -\* b )(y b \\\<8. (5.9) 

Since (recall E(oo) is the ground state energy of H b ) is a Riesz projection of H b we 

have that 

\* b )(* b \ = ^-£(H b -z)- 1 dz, (5.10) 

where T is a circle around the ground state energy of H b in the complex plane with radius 
equal to half of the distance between the first excited state energy and the ground state 
energy of H b . Therefore, 

(\V b )(* b \) S = ^- £((H b - zy l ) 5 dz = ±- £(H bS - z)-\ (5.11) 

where H b5 = (H b ) s . 

Now, we estimate the difference of the resolvents in (|5.10p and (15. lip . The relations (15. ID 
and ()5.3p imply that 

WiHb-HbsXHb-z)- 1 ]] <5. (5.12) 
This and the second resolvent formula 

(H bS - zy 1 = (H b - z)" 1 + (H bS - z)-\H b - H b5 )(H b - z)~\ (5.13) 

imply that, for 5 sufficiently small, 

||(flw - z)- 1 - (H b - z)" 1 !! < 5\\(H bS - z^l (5.14) 

This together with the triangle inequality gives that for 5 sufficiently small ||(i^65 — £) -1 || < 
2||(-ff{, — 2) ||, which together with the previous relation implies that for 5 small enough 

\\(H bS - zy 1 - (H b - zy l \\ < 5\\(H b - z)-\ (5.15) 

Since any z € T is in the resolvent of H b we have that \\(H b — < 1 which with (|5. 15[) 

()5.10p and ()5.1ip implies (I5.9p , This completes the proof of ()5.8p . (|5.5[) follows immediately 
from fl5?3) and (^8l) . 

To prove (|5.6p we observe that similarly to (I5.7P we have 

imps -p)\\< E \\ H ' (m)9 bc (*c\)8 - i^)/ c (^ c i) ii. 

b,ceA at 

Let d > be a constant such that d + H b > 1 > 0. Inserting (d + H b )~ 1 (d + H b ) into each 
summand and using that \\Hs(d + H b )~ 1 \\ < 1, by (|5.ip and (|5.3p . we obtain that 

\\H 5 (Ps-P)\\< £ \\(d + H b )((\* b )g bc (* c \) 5 -\* b )g bc {* c \)\\. (5.16) 

b,c£A at 
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Applying the Leibnitz rule and using that H b ^ b = E(oo)^> b and the triangle inequality we 
obtain that 

\\(d + H b ) ((|^ 6 ) 5 6c (* c |) 5 - |* 6 )ff 6c (* c ) || 

< || ((d + E(^))(\^ b )g bc ^ c \) s - \y b )g bc (* c ) || 

+ \\\V{e- s ^ x - ya) ) ■ Vm b )g bc {m c \e s ^ x - ya) \\ 

+ \\A(e- s ^ x - y ^)\y b )g bc (y c \e s ^ x - ya) \\. (5.17) 

for all b, c G *4 at . Using that, by (13.391) . 93(2; — y a ) > min{P — 1, (2; — y a }}, and using the 
triangle inequality for the bracket (•) and (y a — y b ) > R, we obtain that 

y(z - y a ) + (cc - yb) > R - 1, 

which, together with the exponential bounds (j2.9|) on Vf j, n = 0,1,2 implies that, for 
5 < 9/2, 

|| (x - y a ) rn e- 5v{x - ya) V n ^ b \\ < e" 15 ^- 1 ) , n = 1, 2. 
Furthermore, since </?(x — y a ) < P, we have ||\f r c e l5¥ '( :c_?/a ) || < e^. The last two relations give 

\\\V n (e- 5v{x - ya) )-V m y b )g bc (ty c \e Sv(x - ya) \\<5, n > 0, < m < 2. (5.18) 

The latter equation together with (|5.8p . ([5]) and (I5.18j) implies 

||(d + P 6 ) ((|* 6 }/ c (* c |) 5 - |vP 6 )<7 6c <*c) HI < 5, (5.19) 

for all b, c E A at . This together with (pTTBj) gives ((SSI). □ 

Now, we prove (|3.4ip . i.e. \\(Hg — P)^ 1 ! < 1. We use the decomposition P^- — E = 
H ± -E+ (Hg- - H L ) and that, by {335}, we have that 

IK^-E)- 1 !! <i. (5.20) 

Pulling (.ff — E) out, we obtain that 

Hj--E = (I+ (Hi~ - F ± )(F ± - E)- 1 )^ - E). (5.21) 

We show that (J + (P.5 — H^^H 1 - — E)~ l ) is invertible and estimate its inverse. First, since 
P- 1 := P ± HP ± and := (P" 1 ^ and (P/ - P" 1 ) = P - Jfr, we obtain that 

Hi -H ± = -P 5 ± H S (P 5 - P) + Pt{H& - H)P L - (Pa - P)PP X (5.22) 

Next, we observe that by (|5.ip and (|5.3p . 

||(P 5 - P)(P - P)- 1 )! < ||(P 5 - H)(l - A)- 1 )) < 6. (5.23) 

This, together with flEE) and ([53]) and ||PP ± (1-A)- 1 || < 1 and ||(1- A)P ± (1- A)" 1 ]] < 1, 
gives that 

||(P 5 ± -P ± )(1-A)- 1 || <£, (5.24) 

which for 5 small enough yields \\(Hj- - P ± )(P ± - P) _1 || < 5. This shows that (I + (P^ - 
P- L )(P- L - P)" 1 ) is invertible and its inverse is bounded by 2, which together with (|5.20p 
gives, for 5 small enough, the relation (|3.4ip . 
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6 Proof of (GLE3) and 



Recall that E{oo) = YJj=i E o is the 

sum of the ground states of the atoms. The definition 
(|3,54p of a^p , the relation H = H a + I a and the second resolvent formula, 

(H 1 - - E)~ l = - ^(oo))" 1 - (2^ - E)-\li -E + E(oo))(H^ - E(oo))-\ (6.1) 
where '■= P^H a P^, give that 

<0 = D- G, (6.2) 

where 

D := (fA, (H± - E(po))- x fA)- (6-3) 

and 

G := (f a V a , (H 1 - - E)~ l {li -E + E{oo)){Hi - i^oo))- 1 /^}- (6.4) 
We will show first that 

\G\ < ^ (6.5) 

By (pL26D we have that |E - £(oo)| < which together with (|671j) and the fact that P ± 
commutes with (H — E)~ l and {H^ — ^(oo))" 1 yields that 

G = (P ± {f a ^a), {H L - E)~ x I a {H^ - E{ (X) )y 1 P ± {f^ a )) + O(ig), (6.6) 

Now, by (137101) and {22}, we have that P{f a ^ a ) = # a (*a,/a*a)- But since by (13321) / Q for 
say a = (ran) is odd in the variables z sm := x s — y m for any s £ A m and l^a) 2 is even in the 
same variables we obtain that (faj/a^a) = and therefore P{f a ^a) = 0, which, together 
with P 1 - = 1 — P, gives in turn 

PHfM = /«*«, (6.7) 

for any pair a. The relations (|6.6p and (|6.7p and the boundedness of I a {H^ — ^(oo)) _1 yield 
that 

G = (f a V a , (H 1 - - E^UHi - £(oc))-7/3*a) + O(ig). (6.8) 

Inserting e <5v3 ( x-?/tl )e - ' 5 '^ x-J/a '', where recall that was defined in (|3.39p . into two places 
in JEHU and using that e -M*-J/a) (#1- _ E yl e Sf{x-y a ) = ( H ±. _ where recall #± . 

e -8<p(x-y a ) H ± e s v (x-y a ) ^ see ^jgft), we obtain that 



G = (eM'-vJfA, {Hi - E)- l x f^ a ) + O(^), (6.9) 

where 

X = e- s ^-v«)l a {Hi - E{^))-\ (6.10) 
By ([IS]) we have for S < § that 

|| e M*-!/«)/ a * a || < 1. (6.11) 
By (|6.10p and the triangle inequality we have that 

llxll < \\e- S ^ x - ya \^ ya) <fa(Hi-E{oc))-^ 
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\\e- S ^-^ X{x _ ya} > § I a {H£ - E(oo))^\\. (6.12) 

When (x — y a ) < ^ the electrons are close to the corresponding nuclei and therefore since I a 
has interaction terms between different atoms we have that \I a \ < Therefore, 

1 



e -S^-y a ) X{x _ Ia{H ± _ E^yi^ < J_ (6 . 13) 



On the other hand by boundedness of I a we have that \\I a (H^ — E(oo)) 1 || < 1 and since 

II —Sip(x—y a ) II ^ —cR 

(because ip(x — y a ) > -j when (x — y a ) > -j), we obtain that 



- 5 « {x - Va) X ix - Vn)> nI a {Hi - E{^))^\\ < e~ cR . (6.14) 



-(x-ya)> 



Relations (I6TT21) . ([6TT3]) and (i6TT4"|) imply ||x|| < O(^). which, together with (J6J)]), (l6TTT|) and 
d33lD , yields fT5j) . 

Now, let P a be the orthogonal projection onto *$> a . We show 

D = (/ a * a , {P^H^ " ^(oo))" 1 //?^). (6.15) 

Using the second resolvent formula we obtain that D = (f a ^a, {P^HaP^t — -E , (oo))~ 1 //3\E'a) + 
E, where 

E := (/ a * , " E(^))~ l (P^ H a P^ - P ± H a P ± )(H^ - P(oc))" 1 /^ a ). 

Using that Pj 1 commutes with # a and that H a P ± (H^ - E (oo))- 1 f^a = H a (H^ - E (oo))' 1 
(following from the fact that P 1 - commutes with (H^ — P(oo)) -1 and from (|6.7p ) and 
Pj- - P 1 - = P - P a , we obtain that 

E = (f a *a, (PaHaPa ~ ^M^P - P a )H a (H^ - P(oo))-V^a). 

Inserting e &lplyX ~ ya ^ e~^ x ~ Va ^ twice in the last relation we obtain that 

E = (e 5 ^-y^f a y a ,e- 5 ^ x -y«\P^H a P^ - P(oc))- V^-^ Xl >, (6.16) 

where 

XI = e-M*-*\P - P a )H a (H^ - E{oo)r l f^ a . (6.17) 
Now we will estimate the terms on the right hand side of (|6.16p . By (|2.9p we have 

\\M x -vJf a * a \\ < 1. (6.18) 
Also proceeding similarly to the proof of (|3.4ip . we obtain 

\\ e S>Pix-y a )^p± Ha p± _ J E(oo))- 1 e ^C^-fa)|| < i. (6.19) 

Now we will estimate xi- Since H a = + H a P + PH a + PH a P and the operators H a P, 
PH a , PH a P are bounded we obtain that 

Hfla^-^cx)))- 1 !!^!. (6.20) 



a 
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Also using (|3.10p and (|2.9p . we obtain that 

e -M*-?/a)(p_p a ) = ^ e- 5v{x ~ Va) \^ b }(^ b \ = 0. (6.21) 

The last two relations and (|6.17p give that ||xi|| = 0. This relation and relations (|6.16p . (|6,18p 
and (|6.19p imply E = 0, which gives (|6.15p . 

We will now consider the first term on the r.h.s. of (|6.15p . Consider an atomic de- 
composition a = {Ai, Am} and let a = ij,/3 = kl. Recall that the atomic hamiltonians 
H Ai and eigenfuctions (pAi were defined in (|2.4p and (|2,7p . respectively, and recall that Ei 
denote the ground state energy of the atom Ai, so that (omitting identities) H a = Yli^Ai 
and \l/ a = \\ i 4>Ai (see (|2.3p and (|2.6p ). Let Pi = \4>Ai)(4>Ai\, the orthogonal projection on 
4>Ai > and Pij = Pi (g> Pj and P lJ = Yl^ijPk- Without changing the notation we consider 
the operators P^i and P kl on the entire Hilbert space of the system and, abusing the 
notation, write the orthogonal projection P a on P a as P a = Pki <8> P kl ■ With these notation, 
we have P^~P kl = P^ ® P kl and therefore, on the invariant subspace RanP^, 

i^k,l 

where Hfc = Pjti(HA k + H Al )P^. This and the relation f M ^ a = P kl f M ^>a give 

{PtH a Pt - J B(oo))- 1 /fc^a = {PtliHil ~Ek- E^Piif^a, (6.23) 

which implies 

(/«*«, {PtH a Pt - E{^))~ l f^ a ) = W, (6.24) 

where 

W = (fi^ a , (Pii(HA k + H Al )Pti ~Ek- £?i)-Vfcl*a>. (6-25) 
f if ij ^ kl 

We now prove that W = < ' .„ . _ , . , where are given in (|3.43p . If ij / kl, 

y crjj, it % j — k,i 

then we may assume without loss of generality that i ^ k,l. Observing that the inner product 
on the right hand side of (|6.25p reduces to an integral that is odd in z m i,m € Ai we obtain 
W = for ij ^ kl. For the case ij = kl, W = aij follows immediately from the relations 
(|2.6p - (|2.7p and the fact that the ground states of the atoms are normalized. This, together 
with (jO]) . dSSJ, USE), (523) and (IBT25]) . proves the relations (jSMj) and ((336]). 



7 Proof of Theorem 11.21 

7.1 The general set-up in the case of statistics 

We begin with the general setup for the proof of Theorem 11.21 We will follow closely Section 

CO 

Let Sn denote the permutation group of {1,2, N}. We consider the representation 
T : S N -> U(L 2 (R 3N )), where U(L 2 (R 3N )), denotes the unitary operators on L 2 (R 3N ) defined 
by the relation 

{T^){xi,...,x N ) = ipix^-i^^x^-i^, ...,x w -i {N) ). (7.1) 
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We label by a the types of the irreducible representations of T. Let Ti 17 be the subspace of 
L 2 (]R 3Ar ) on which the representation T is a multiple to the irreducible representation of type 
a. Then 

L 2 {R 3N ) = ® a H a . (7.2) 

Let Xg denote the character of the representation of type a evaluated at g. The orthogonal 
projection of L 2 (]R 3Ar ) onto 7i a is given by 



^ = tf E X^-iT*, (7-3) 

7rs5jv 

(see |Haj ). Since H commutes with any T w , by (|7.3p . so it does with Q a , 

HQ a = Q a H. (7.4) 

Hence T~L a is an invariant space of H . Let H a be the restriction of H onto W . We now 
denote the ground state energy of the system for the irreducible representation a by 

We repeat some definitions from the introduction. Fix an irreducible representation Tg 
of type a of the permutation Sn- Let S(a) be the subgroup of the permutation group Sn 
that keeps the clusters of the decomposition a invariant. Then the restriction Tg \g^ of the 
representation of SV of type a is not necessarily irreducible. Therefore, there exists a family, 
7°", of irreducible representations of S(a) such that 

T S N \s(a) = ®a£l°Tg( a y 

The representations a € I a are called induced representations and we write a -< a. The 
corresponding decomposition a will always be clear from the context. 

We denote by the orthogonal projection onto the subspace of L 2 (M 3Ar ) on which the 
representation of S(a) is a multiple to the irreducible representation of the type a. By [Ha] 
it can be written as 

«-#t)E^ 9 , (7„ 

g€S(a) 

where #S(a) is the cardinality of S(a). Since H a commutes with all permutations in S(a), 
it commutes with Q" . We denote by the restriction of H a onto Ran ■ 
We write a a if a -< a and inf <r(H") = min^^ai inf a(H^). Let 

E a (oo) := min mfcr(H^). (7.6) 

Let be the orthogonal projection onto the ground state eigenspace of and P a be 
the orthogonal projection onto the subspace 

span{RanP"|a G A at ,a -<-< a}. 

We show below that P a commutes with T n , for any ir € Sn, and therefore with Q a and 
Ran P ff n RanQ 17 ^ {0}. Since P a and Q a commute, the operator 



IT := P a Q a = Q a P a , 



(7.7) 
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is the orthogonal projection on Rani 30 " D RanQ CT 7^ and satisfies 

WQ U = Q U W = IF. (7.8) 



We now consider the Feshbach map, Fr<t(X), where 11°" is defined in (|7.7p . as defined in (|3.1 
- D, i.e. 



F n *(X) = [WHIT - WHW L (W^HW^ - \)~ 1 W J - L HIT 7 ] | Ran n-, (7.9) 
with, as above, n CT_L = 1 — IT 7 . Since H and i-" 7 commute with Q a , we have that 

F n .(A) = [P°H°P° - U°(\)] | RannCT , (7.10) 

where 

U a (\) = P a H a P a± {H a± - \y 1 P a± H a P a , (7.11) 
with P a± = l-P a and H a± = P a± H a P a± . As before (see (£0])) we have 

A eigenvalue of H a A eigenvalue of Fn^(A) (7-12) 

and to prove Theorem 11.21 we have to estimate the different terms on the r.h.s. of ()7.10p . 
We conclude this section with some general results used above and below. 

Lemma 7.1. P a commutes with T n , for any ir E Sn, and therefore with Q° . 

Proof. To show that T n commutes with P a , it suffices to show that if ^ G Rani 30 ", then 
T w ^f G Rani-" 7 . To show the latter it suffices to take ^ satisfying 

H a m = E u {oo)^. 

Therefore, (T 7r ff a T 7r -i)T 7r 1 I' = E a (oo)T n *$> . In addition we have that 

TnHaT^—i — Hna> 

where the clusters of b = ira were defined in (12. 2ft . The last two relations imply that 

H m T^ = S CT (oo)r^. 

As a consequence, G RanP°". □ 

Lemma 7.2. For each a -<-< a, any normalized <1> G RanP" and for R is large enough, we 
have the estimate 

\\Q a n 2 = (7.i3) 

and therefore Ran P a D Ran Q a ^ {0}. 

Proof. Using that, due to the exponential decay of the ground states of the atoms, $ and 
T n Q, for 7r G S]y/S(a), are almost orthogonal i.e. (§,T n §) = 0,V7r G Sn/S(o), and using 
the expression for Q a , we have that 



= f E x:-i($,T^}. (7.14) 

' 7re5(a) 
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Since P^Js^) = ©a^o-P^), taking the trace of T^_ 1 = ©a^o-P"_i, V-zr G S(a), we obtain 

4-i = ExJ-i 5 V7re5(a). (7.15) 



This, together with the expression (|7,5p for the projection and the relation (|7.14p gives 

that 

The last relation together with the fact that Q^^a' 3 = bafi^a' 3 gives (|7.13j) . □ 

For given a € A at and a -<-< <r, let 

{$^',j = l,...,deg(a)}, (7.16) 

be an orthonormal basis of RanP", where recall that P" is the orthogonal projection onto 
the ground state eigenspace of H®. Maybe we should mention it here. Here deg(a) is the 
degree or dimension of the representation a (recall that from Condition (D) the ground state 
eigenspace P" of H% consists only of one copy of the irreducible representation of type a). 
By the last lemma we use T^a' 3 = &na ■ 

Lemma 7.3. Each basis functions $a'* are products of ground states of atoms: 

= (7.17) 

where i are identified with (i\, . . . , im)- 

Proof. We have that S(a) = (^jL 1 S(Aj), where Aj are the clusters of the decomposition a, 
and S(Aj) is the permutation group of the set Aj. Since a -< a, a is of the form a = ®Y=i a j 
where ay is an irreducible representation of S(Aj). To show (|7.17p it suffices to show that 

AI 
3=1 

where recall that P^ denotes the projection onto the multiple of irreducible representations 
of type p. We have that S(a) B n = m...TTM, with 7Tj E S(Ai), and 

M 

#S{a) = J] #S{Aj) and T n = T ni ...T nM . 

The last relation and the definition of characters imply that 

M 

The last two relations, the formula (|7.5p and the relation x°-i = X°-i & ve (|7.18p . □ 
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Lemma 7.4. The one-electron density, 

dcg(o) 

pi(*i):= \^' l (z 1 ,...,z N )\ 2 dz 2 ...dz N , (7.20) 

i=i J 

written in terms of the variables z defined in (|3.19|) . is spherically symmetric. 
Proof. For any rotation R we define the action 

T R <f>(z 1 ,...,z N ) = <!>{R- 1 z 1 ,...,R- 1 z N ). 
Then a simple change of variables together with (I7.20|) gives that 

deg(ct) 

\T R <PZ i (z 1 ,...,z N )\ 2 dz 2 ...dz N =p 1 (R- 1 z 1 ). (7.21) 

i=i ^ 

Also since the Hamiltonian iJ° commutes with Tr we have that T R & a ' 1 is also a ground state 
of H£ and therefore by Condition (D), 

deg(a) 

T R ^ i {z 1 ,...,z N )= d ik ^ k {z u ...,z N ), (7.22) 

k=l 

where (da-) is a unitary matrix. Therefore, by (|7.2ip and (|7.22p we have that 

deg(a) deg(a) 

p 1 (R~ 1 z 1 )= Y \ Yl d ik ^ k { Zl ,...,z N )\ 2 dz 2 ...dz N . (7.23) 

i=l J k=l 

Expanding the right hand side of the last equation and using that the matrix dij is unitary 
and therefore satisfies X^Sl dik 1 dik 2 = <5fcifc 2 ; we obtain that 

dcg(a) 

p 1 (R- 1 z 1 )= Y m' k (z 1 ,...,z N )\ 2 dz 2 ...dz N . (7.24) 
k=l 

The relations (|7.20j) and (|7.24j) imply the spherical symmetry of pi(z±). □ 



7.2 Estimate of P a H a P a 
Lemma 7.5. 

P a H a P a = E a (oo)P a Q a . (7.25) 

Proof. In this proof we omit the superindex a in P a , E a {oo), Q a and write instead P, E(oo), Q. 
Recall that F° is the orthogonal projection onto the ground state eigenspace of H®. Using 
the orthonormal basis f|T. 16[) in RanP", we write 

deg(a) 

P:= E (7-26) 
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Due to the exponential decay of ground states of atoms together with the fact that Ran P", Ran P' c 
are orthogonal whenever a ^ j3, Va,/3 ~< a, and the fact that &a' J ,j = 1, deg(a), is or- 
thonormal we have that 

($<M = SabSapSij, Va, b G „4 af , a, /3 ^ a, i = 1, . . . , deg(a), j = 1, . . . , deg(/3). (7.27) 

Using the last relation together with (|7.26p , and the fact that P is by definition the orthogonal 
projection onto span{RanP"|a G A at ,a -<-< a}, we obtain that 

deg(/3) 

p = E p "= E E i*fx*? J i- ( 7 - 28 ) 

&G.A at ,/3^(T 6e^l ai ,/3^^(T i=l 

Using (17371) . (I7T28I) and the relations H = H b + I b and # 6 $f j = J3(oo)$£ J ', V6 € 
A at ,/3 -<-< er, we obtain that 

PfTP = E{oo)PQ + UQ, (7.29) 

where 

E I^X^'S^'X^'I. (7-30) 

a,b,a,/3,i,j 

We will now show that 

V = 0. (7.31) 
If a ^ b, then <3?a'\ have exponentially small overlapping tails and therefore 

(*«-', J 6 ^') = 0, Va, 6 G A at , a ^ 6. (7.32) 

Let a = b but a / Since I a commutes with for all permutations in 7r G 'S'(o) and 
therefore, due to (I7.5p . with Qa, we have that I a &a € RanQa- Since also G RanQ", 
we obtain that 

lo^) = 0, Va, /3 -<-< a with a^. (7.33) 

Assume that a = b, a = (3. Let c 71 " be the matrix of TV, ir G 5(a), w.r.t the to the basis (|7.16p of 
Ran P". Then, since T w is unitary and commutes with I a and due to 1^$"'* = Ylt^i^ c ik^a' k 
and similarly, T n ^' j = Zflf^ , we obtain that 

deg(a) deg(a) 

= E E ^(^W). 

fc=l Z=l 

So if we define the matrix h with hji = (&a' 1 , I a &a^), then the last equality becomes 

deg(a) dcg(a) 

h i i= E E °ji h ikWk, 

s=l t=l 

or h = c^hc 77 *. Since = 1, we obtain that hc n = c n h, Mir G S(a). Since the repre- 

sentation a is irreducible we have, by a standard result of representation theory, that h is 
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a multiple of the identity. (Indeed, let u be an eigenfunction of h with eigenvalue A. Then 
hc n u = c n hu = Xc n u. So the eigenspace of the eigenvalue A is invariant under the represen- 
tation and since the representation is irreducible the eigenspace has to be the entire space. 
So h is a multiple of the identity.) This shows that 

{*? i Ja*F)=8 ij {*Z\l a *? 1 ). (7.34) 
Now, using Lemma 17.41 and proceeding as in the proof of (|3.16p for a = b we obtain 

deg(a) 

£ (^Ja^)=0, (7.35) 

3=1 

which by (|7.34p implies 

Ja*?*) = 0, Va «a,i= l,...,deg(a). (7.36) 

The relations fPOj) . (fT32]l . (fT33]) . (17T341) and (I7T36D . imply (TOTD . From (|7T29j) and <|7^I|> 
we obtain ([7725]) . □ 



7.3 Lower bound on H aL 

Let E a (oo)' denote the first excited state energy of the system of non interacting atoms and 
let 

7 f= mm (infa(F a Q )-^(oo)), 72 = E' T (oo)' — E a (oo). (7.37) 

a€A/A at ,a^a 

By Property (E) we have that 7J" > and, by the definition, we also have that 7J > 0. We 
also define 

7o CT = min{ 7l \7 2 T }- (7.38) 
Lemma 7.6. For R sufficiently large, there exists 7 CT > such that 

H a± > (E a (oo) + 2 7 <T )Q <T . (7.39) 

Proof. In this proof we will also omit for simplicity the superindex a in P a , Q a , E a (oo), 7J 
and write, instead, P,Q,E(oo),jj. Using that H 1 - = P- L H cr P- L and that Q commutes with 
P 1 - and the relation (|4.9p and P 1 - < 1, we obtain 

P ± H ap± > QP ± J a H a J a P ± Q - ^Q. (7.40) 



R 



We will now show that 



QP ± J a H a J a P ± Q > (P(oc) + 7o )QP ± JfP ± Q - ^Q, Va £ A (7.41) 

it 

Case 1: a € ^4/^4 a *. By (17371) and f7T38l) we have that 

^ = P?H a > (E(oo) + 7o )Pf (7.42) 
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Using this, the relation Q = QYlp^Qa and the fact that Q a commutes with J a and Q- 1 , 
we obtain that 

QP^J a H a J a P^Q > (E(oo) + J )QP ± J a E P^aP^Q 

= (E(oo)+ l0 )QP ± J 2 a P ± Q, 

which gives (|7.4ip for a G A/A at . 

Case 2: a G A at . Recall that P" is the projection onto the ground state eigenspace of 
H£ and define 

Pa := ®a^aPa- (7-43) 

If b € A at with b ^ a, then, repeating the arguments of the proof of (]4.15p . we obtain 

\\^ j J a H a \\, \\H a J a ¥ b ' j \\ = 0, Vj = 1, ...,deg(/3). (7.44) 

Relations (THZBD . (17^31 and (fT4ljl imply that 

H a JaP b = 0, P b J a H a = 0, V6 € .4°' with 6 ^ a, (7.45) 

and therefore, since P = J2beA at anc ^ ^ >_L = ^ ~ P> -^a" = 1 ~~ Pai we obtain 

QP ± J a H a J a P ± Q = QP^J a H a J a P^Q. (7.46) 
In the same way as (|4.22p . we show that 

||[|^)(^'|, J a ]H a \\ < 1 Va ^ a, Vj = l,...,deg(a) (7.47) 

and 

||F a [J a ,|^)(^|]||<i, Va^a, Vj = l,...,deg(a). (7.48) 
By (T7T26]1 and (17151) we have 

deg(a) 

P a= E E I^Wl, 

a^Xcr j = l 

which, together with P^ = 1 - P a , fTTTP and (fTTigj) gives 

P a L J a H a = J a P^H a + 0(^) and H a J a P^ = H a P^J a + O(^). (7.49) 

From (|7.49|) and the fact that Q = Q Ylp^a Q°- an< ^ commutes with Pj 1 and J a for all 
/? -< a, we derive 

QP^JaHaJaPtQ > QJ a Pa E H%P£j a Q ~ ^Q. (7.50) 



If /3 -< cr but /3 7^ a, then by fTBTP we have that 

fl£ > (P(oo) + 72 )Qf , /3 -< a, (3 M (7.51) 
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If now j3 a, then by (|7.43p and = 1 — P a we obtain that P^Qa projects out the 
ground state eigenspace Ran Pa of Ha and therefore by ([7.37P , 

PtH^ > (E(oo) + J2)P^Q% Pa, P ^ (7-52) 

From relations (|7.38|) . (|7.5ip and (|7.52p we obtain that 

P a L H?P a L > (E(oo) + 70 )P a ± Qf P^, V/3 -< a. (7.53) 

Repeating the arguments of the proof of (|7.49p we can obtain that 

In addition, we have that 

JaP a Ja = P a Ja^a + ['^aiP a ]^a ~ P a Ja\JaiP a ]• 

The last two relations give us that 

C 

j p ± j < p 1 - j 2 p ± -i 

' J a r a ,J a 2i r a J a r a ' p' 

Using the last relation and the fact that 

P^Ja = P ± Ja, JaPa = ^ , Va € A*, 

which can be proven similarly to (|7.45|) . we obtain 

JaP a L Ja<P ± J 2 a P ± + ^- (7.54) 

Similarly to (|4.28p . we can show 

E(oo) + 70 < 0. (7.55) 

The relation ([711")) for a G A at follows immediately from the relations ([Tiff]) . ([73U]) . ([733"]) . 
(|7.54p and (j7.55p . together with the relations Q = Q J2p <cr Qa and Ha = QaH a and the fact 

that Qa commutes with P^ and J a , for all /3 -< a. This completes the proof of (]7.4ip . The 
relations gZD, ([710]) . ([71T]) together with (1735]) and Q 1 - < 1 give that H 1 - > (E(oo)+^ )Q, 
provided R is sufficiently large, as desired. □ 



7.4 Rough bound on E a (y) 
Lemma 7.7. For R sufficiently large, 

E a (y) < E a (oo) + 0{e~ CR ). (7.56) 
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Proof. We fix an a a and let := ^a' 1 . Since Q° is self-adjoint and, by (|7.4p . commutes 
with and since, by Lemma 17, 21 ||<2 CT< £a|| > 0, we have that 

Therefore, using the last relation, that H = H a + I a and that H a ^" = £' CT (oo)$°, we obtain 
that 

E°{y) < E°{oo) + W * a l a l . 

\\Q a ^a\\ 

It is therefore enough to show that 

~Wn\T ~ ' ( } 

Due to exponential decay of ground states we have that 

(T^S, I a $«>=0, Vn€S N /S(a). 
With this observation, we can repeat the arguments of the proof of (I7.13P to obtain 

The last relation together with ([77l3|> gives that ^qJ^I"^^ = ($£>4$£) and therefore by 
(I7.36P we obtain that 

The last relation gives (|7.57p immediately. □ 



7.5 Estimate of U u := 

We show in Appendix IA1 that the term U a (E a ), defined in (|7.11|) . satisfies the estimates 

1,M o,a 1,M 



where 



cH-«r i<j W i<j W 



P a := Q CT ^ P«Q° (7.59) 



a£A at 



and <r^' a are given in (|A.4j) . 
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7.6 Completion of the proof of Theorem 11.21 

Since IF = Q tJ PQ (J by (TT28]) and (f739l) we obtain that 

if = pa - ( 7 - 60 ) 

Therefore, from relations (I7TT01) . ([735t (|735|) and fTpO]) we obtain that 

(1,M a,a l.M \ 

Moreover, using an argument similar to the one of the proof of the relation (|3.36p we obtain 
that E a is the lowest eigenvalue of F n * (E u ). By (fT59|) and by the fact that P^pff = 0, Va ^ P 
and that commutes with Q a for all a € A at ,a a we obtain that P a P@ = for all 
a ^ (3 which together with (|7.6ip and the fact that E a is the lowest eigenvalue of Fu<t(E' 7 ) 
gives that 

1,M cr,a 1,M 

E° = E°(co) + min {- £ — ^ — ^ + 0(£ 

as desired. 



A Estimate of U a := U a (E a ) 

Recall that 0^' fc ,j = 1, n m , form an orthonormal basis of the ground state energy space 
of the m— th atom with the ground state energy E° n , 

and that as in ([22]), (t> a A .{x Aj ) = ^' % (x A . with XA i = ( x i '■ i ^ x Ak -yj = (xi~yj ■ 

i E Aj) and (j)^' 1 a ground state of the j-th atom with the nucleus fixed at the origin and of 
a symmetry type a. Let P± := 1 — Pij and Pij = Pi® Pj, with 

and 

^ := + H Al )Pii - El - E[)-\ (A.3) 

Our goal is to prove the following lemma: 

Lemma A.l. The equation (I7,58p holds with 



■Hi TLj 



'= dik) E E (fi^AM^^(H Ai + H Aj )P± - El - Ejr^f^l). (A.4) 

' l=\ m=l 



Long Range Behaviour of van der Waals force, April 22, 2012 



42 



Proof. In this proof we will also omit for simplicity the superindex o in i* 7 , Q a , E a (oo) , E a (oo) , j a 
and write, instead, P, Q, E, E(oo), 7. Recall that U a (\) was defined in (17.111) . Using that 
H a = HQ and the fact that Q commutes with H and P we obtain that 

U a = QPHP- L (P ± H a P ± - Ey 1 P ± HPQ. (A.5) 

We extend the class of decompositions by considering, in addition, partitions into a molecular 
ion and the j— th electron and define the refined covering 

Of := {\ Xl - y k \ < f3R,Vi G A k ,Vk = 1,...,M}, Va G A, 

:= {x G R 3N : \x, - y k \ > pR/2,Vk = 1,...,M}, Vj = l,...,N, (A.6) 
and the corresponding partition of unity { J a , Jj} having the properties 

< J a < l,suppj a C nI,||VJ a ||L« < -^,\\AJ a \\ L ~ < Va G AU{1,2,...,N}, (A.7) 
and obeying the relation 

ae^U{l,...,JV} 

The construction can be proven following |Sig| as we did for the construction of the partition 
of unity satisfying (|4.5p and (|4.7p . As in Subsection 17.31 we represent the operator P- L H cr P- L 
as 

P ± H CT P ± = QKQ, (A.9) 

where 

K = K + O(^), (A.10) 
6e.4U{l,...,iV} 

K . = J b PtHgP b ^J b , for 6 G A 

[JjEp^Hfjj, torb = j€{l,...,N}, 

where if 1 " = 1 if 6 G ^4/^ at or /3 a and if - 1 = 1 - if if b € -4 a * and /? -<-< a, 
and iij := H + X^fcf=i \x--y k \ ' Moreover, for 6 = j G {1, A r }, /3 refers to an irreducible 
representation of Sn-i- 

Lemma A. 2. There exists 7 > suc/i that 

K>E + 2j, K >E + 2j. (A. 13) 

Proof. Due to (|A.12p and the HVZ theorem with symmetries (see |Sig| ) we have that there 
exists 71 > such that 

Kj >(E + ji)Jj, Vj G {1, N}. (A.14) 
In addition, by ()7.53p and ()A.12p we obtain that 

K h > (E(oc) + 70) E J b P b L Ql P b Lj b, V6 G A. 

8<a 



with 
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But since E(oo) + 70 < and 



we obtain that 

K b > (E(oo) + 70) Jl V6 G A. (A. 15) 

Therefore, by (|A~8|) . dATTOl) . ([Axil) . (jA~T4|) and (jA~T5|) we obtain (jATT3|) . □ 

From ()A.5P and (|A.9|) we obtain that 

C/ CT = QPHP ± (QKQ - E)~ 1 P ± HPQ. 

Since K is symmetric with respect to particle coordinates, it commutes with all permutations 
and therefore with Q. Since Q commutes with P, H and K and due to (|A.13p . we have that 

U a = QPHP ± (K - E)~ 1 P ± HPQ. 

On the other by the second resolvent formula we have that 

(K - E)~ l = (K - E)~ l -(K- Ey l (K - K )(K - E)~\ 

The last two relations together with (|A.10j) imply that 

U a = QPHP ± (K - E)~ 1 P ± HPQ + Rem, 

where 

Rem = QPHP ± 0(^-)P ± HPQ. 
Therefore, using the last relation and 

p= E i^'x^'i, 

adA at ,a^,^a,j=l,...,dcg(a) 



we obtain that 



deg(o) deg(^) 

U<T = E E E Q^VTH^IQ* + Rem, (A.16) 

a,b&A at ,a,fi«cr i=l j = l 



yaPij = ( H <S>%\P L (K - Ey 1 P ± H^ b ' j ). (A.17) 



where 

Using now that H = H a + I a and that is eigenfunction of H a we obtain that \HP 
(I a $a' j \P- Similarly we obtain that P^H^ = P ± I b $ b 3 ' j . This gives that 



V°f ij = (I a ^ i ,P ± (K - Ey 1 P ± I b ^). (A.18) 
Repeating the arguments of the proof of (13.4T[) . we show that 

vT j = <W, (K - E)~ l I b ^). (A.19) 
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Similarly we can obtain that 

deg(a) deg(/3) 

Rem= E E WMr)^)- 

a,beA at ,a,/3^-<cr i=l j=l 

Proceeding as in the proof of (|3.30p we can obtain that 



\\Ic^ k \\ < 0(E | _! | 3 )>Vc € A at , 7 ^ a,k = 1, ...,deg( 7 ). (A.20) 
The last two relations give that 



Rem = 0(^-_L- i7 ). (A.21) 



The relations (|A. 16|) and (|A.2ip together with Lemma IA.3I below and (|7.26p . gives (|7.58p ~ 

Lemma A. 3. T/ie following equations hold: 

1,M 

r a(3ij 



V^i = 0) Va G .4 at , Va, /3 ^ a, a / /3. (A.23) 

Ca^ = "C all %- (A.24) 
F a °« n = VJg* 11 , Va, b G ^ a *. (A.25) 



Moreover, we have that 

1,M 4 cr,o 1,M 

tzj\y*-yi\ e ■"■ 

where the positive constants a^- a are given by ()A.4p . 

Proof. In this proof we keep writing E and E(oo) for and E a (oo). For each c € A at we 
consider a C°° function Xc that commutes with all permutations in S(c) such that 

*,(*) = { lifXG 5 (A-27) 
{ if x € (ft c 2 ) c , 

where J7| was defined in (|A"|) . and 

l|9 Q X c || < #~ H , for any multi-index a. (A.28) 
Then due to the exponential decay of <3?c J we have that 

4$™ = xdc^c\ <&™ = xM J . (A.29) 
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The last relation and ()A.19P imply that 

Vf ij = (I a ^\ Xa (Ko - Ey^I^). (A.30) 

Since XaXb = and 

[Xa, (K - E)- 1 } = (K - E)- l [K , X a]{Ko ~ E)~\ 

we obtain that 

Xa(K - E)~ 1 Xb = (Ko ~ E)~ x [Kq, Xa]{K - E)~ l Xb . 
Furthermore, since XaJc = 0,Vc ^ a the relations (lA.lip and (|A.12p give that 

[Ko,Xa] = [Ka,Xa\- 

The last two relations give that 

Xa (K - E)~ 1 Xb = (K - E)- x [K a , Xa](K - E)' 1 . (A.31) 

Also by (jA.lip we obtain that Kq is symmetric with respect to the particle coordinates and 
therefore 

T 7T K = K T n ,WeS N . (A.32) 

Hence, due to (|7.5p . Kq commutes with Q". Since Q"/ a = I a ^a' 1 and commutes with 
K a ,J a and K , the equation $A~3l} together with (jATT2|) . (fAT30|) and Q%Q p a = 5 a ^Q%, gives 
that 

vT j = [J a H2P^Ja,Xam, (A.33) 

where 

* a = (Ko - Ey 1 ^ 1 , V& 6 = (K - E)- 1 !^. (A.34) 

Since H^P^ 1 - = H® — E^P" (where E% is the eigenvalue corresponding to the eigenspace 
P° of H") and Xa commutes with J a we obtain that 

[J a H2P: L J a ,Xa] = Ja([H«,Xa] ~ E«[P« , X a]) Ja- 

Due to (|7.26p and ()A.29P we obtain that [P" , Xa] = and therefore 

[JaHaP^ ± J a ,Xa] = J aWa i Xa]J a ■ 

Using now that H% = H a Q" and that Q% commutes with Xa and using that [H a , Xa] = 
[-A, Xa] = -Axa - (Vxa) ■ V, we obtain that 

[JaHaPa L Ja,Xa] = J a(Qa*[Pf a , Xa]QZ) J a 

= -J a Q«(A Xa + (Vxo) • ^)Qa Ja- (A.35) 

Using the last relation, together with (fA"l3|l . (1AT20]) . (1AT28]) . (|AT33|) and (EOIj) . we obtain 
that 

l,M 

' = -<*., 4Qa((VXa) • V)Q a a )Ja* b ) + 0£ r — , 7 )- 

r ^ ■ \Vi Vj 
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To estimate the first term on the r.h.s., we note that there exists a positive constant c such 
that 



X 1 := HVQ^JAII < \\^/KT~cQ2J a ^a\\ = (^a,JaQ2(H2 + c)Q2Ja^a} 1/2 . 

Since H% = H%P* + H° P a a± and H^P" is bounded operator, by (TOOjl . this gives 

1,M 

X l ^ °(E I 16) + UHZK 1 - + c)Ja*a)- 

Since by (jA.lip and (jA.12p J a H^P^ J a is summand of Kq and each of the summand of Kq 
in (jA.lip is bounded from below we obtain that there exists a constant d > such that 

JaiKP^ + C)j a <K +d. 

The last two relations give that 

1,M 



The last inequality together with ()A.20p and ()A.34p gives 



i,m 

:= ||VQ£J a * a || < 0(]T I is)- (A.36) 

Therefore, <|Xl2)> holds true. Using (lA~l7|) . that Q%$l' m = 5 aj $2' m , and that is self- 
adjoint and commutes with I a ,Q, H, Kq (with the last because of ()A.32p ). we obtain ()A.23p . 
Using ()A.17p and repeating the arguments of the proof of (|7.34p we prove (|A.24p . 

We will now prove (|A.25p . Since a, b S A at there exists a permutation ir such that b = ira 
where the action of n was defined in f)2.2[) . Hence we may assume without loss of generality 
that T^^ 3 = Since on the other hand T v is unitary and commutes with Q a , P, H, Kq 

and T n I a = I^T^ we obtain using (|A.18|) that 

V a a a ajJ = (I b T^y,P ± (K - E°y x P x I b T^). 

But since b = vra and T^f j = $™ ', (|A~25|) follows. 

We will now prove (|A.26p . We first prove the following relation: 

1,M 

Yi := V™" - (I a $y, (K a - E)' 1 ^) = o£ = (A.37) 

j " : -'0 

Using (IA.19P and (jA.29j) . we obtain that 

Y 1 = (I a ^ j , [{Kq - E^Xa ~ Xa(K a - Ey 1 }^). (A.38) 

By ((SI), (IA.7P and (IA.27P we obtain that XaJb = for all a / 6 which together with (jA.lip 
and (jA.12p gives that 
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The last relation together with (|A.38|) and the second resolvent equation gives that 

Y x = (I a $y, (K - E)- l [ X a,K a ](K a - E)- l I a $y). 

Since Q^I a ^2' J = I a ^a' J and Q% commutes with both K a ,Ko the last relation together with 
([SZE1D and Q a a Qi = 0, for j3, gives that 

Yi = {I**? , (K - Ey l [ Xa , JaKP^J a ]{K a - E)- x I a ^). (A.39) 

The relations (jOoll . (EP8j) . (jSMJ) and fOojl yield that 

Fi = (I a *^', (Ao - E)" 1 J a {Q a a {VXa) ■ VQ%)J a {K a - E)~ 1 I a ^^) 

1,M 

Therefore, by (|A.13[) . rtA.20h . (1A.28 j) and the last relation, showing (|A.37|) reduces to showing 
that 

1,M 

\\VQ«J a (K a - S)- 1 W|| < OQ2 | jg), (A.40) 

which can be proven similarly to (|A.36p . This completes the proof of (|A.37|) . 
We will now prove that 

F 2 := (I a ^' j , (K a ~ E)- l I a <$>y) - {I^y, {K« - i?)- 1 W) 

1,M 

where 

K- = P^H-P^. (A.41) 
Using the second resolvent formula we have that 

(K a - Ey 1 = (JC - EY 1 - (iC - E°y\K a - K«){K a - E)~\ 

and the relation Q"/ a $a J = Ia^a' 3 and using that Q" commutes with K%, we obtain that 

F 2 = {I^y, (K« - E)~ l Q%(K a - K2)(K a - E)' 1 ^). (A.42) 

In addition, using (|A.12|) . ()A.4ip and that Q^Qa = 0, we obtain that 

Qa^a = Qa^a,K a J a - 

The last relation together with ()A.42p gives that 

f 2 = {i^y, (k« - E°)-\j a K%j a - K-)(K a - EyH^y). 

Inserting —K^J^ + K^Ja we obtain that 

F 2 = F 21 - F 22 , (A.43) 
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with 

Y 21 = (I a $y ,{K* - E)- x [J a ,KZ\J a {K a - EY l I a $>y) (A.44) 

and 

Y22 = (IaK' j , (K ~ ET'Kil - J 2 a ){K a - E)~ x I a ^). (A.45) 
We will first show that 

x,M 

\Y 21 \<0(£] ir)- (A. 46) 

Indeed, by (TOT]) and the fact that K% = H% - E(oo)P£ we have that 

[J a ,K a *] = [J a ,H a ]-E(oo)[J a ,P a *]. 

Since also [J a , P£\ = (we can obtain that by repeating the arguments of the proof of (|7.49|) ) 
and [J a ,H a ] = [J a , —A] we obtain that 

[J a , K«] = P«AJ a + P«(V J a ) ■ V + 0(1). 

Using (|A.7[) . (|A.20p . (|A.44j) and the last relation, we obtain that showing (|A.46[) reduces to 
showing that 

(I**?, (KS - Er l P2{{VJ a ) ■ V)J a (K - E)' 1 ^) 
1,M 

< v " 

Using (|A.7p and (|A.20p we obtain that the last estimate reduces to showing that 

1,M 

\\VP?(K2 -E^Ia^W < 0(E u,_,„| 3 )- (AA7) 

Since 



\\vp%{K - E)- l iaK' ] \\ < \\VKTc~p:(k« - e)- 1 

we obtain that 

\\VP:(KZ - E^I^f < (I a 1>y, (K: - E°y\H2 + c)(K« - E)' 1 ^). 

The last relation together with the fact that H a = H^P^ + H^P* = K^ + E^P^ <K% + E% 
and the relation ()A.20p gives (|A.47j) and this completes the proof of (|A.46p . 

Now, we estimate ()A.45p . We will now show that Proceeding as in the proof of (|3.4ip we 
can show that 

WdK^s-E^W <1 (A.48) 

(in fact showing the last relation is much simpler than (|3.4ip because P° is eigenspace of 
ff"). Inserting e - s v(^-Va) e 5<p(x-y a ) three times we obtain that 

y 22 = (eM x -*h a $p, {{K2)s ~ E)-\K2)se-M*-»>\l - J 2 a ){K - E)' 1 !^). 
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It follows from from ([339]) . (jAj), (1X71) (supp J a C Ol) and (jOJ) that 

e -M*-»«)(i_ j2) = . 

Since, by (IX48|) . )<5 - E^iK^s is bounded, this leads to y 2 2 = 0. This estimate and 
relations (1X43|) and fOB]) . give ((Aj). The relations (fA~37l) and JA} imply 

i,Af 

Using (|A.49p and proceeding similarly to the proof of (|3.53|) we obtain that 

deg(a) 1,M 1,M deg(a) 1,M 

£ ^— = EEE e4 ^S + o(£ T-^-jf). (a.50) 

m=l i<j k<l m=l i<j ^ % ^ ' 

where 

%rM ' \yi - Vj?\vk - yi\ 3 

By (|7.17p . we have that each <&a' m , ,m = l, ...,deg(a), is product, 

*? m = ^ <1 ...^5£', (A.52) 

of ground states, 4>°^ 1 , of atoms with the ground state energies E% (see (|A.1|) ). Using (|A.52|) . 
we show as in (16.23P that 

((P^H^P^ - Eioo))' 1 ^ = Rbfu***, (A.53) 

where the operators are given by (fO]) . Now, using (|A.53j) . we obtain that 

m = (fv*? m ,Rtifki*sn (A54) 



= 0, when »j / fcl. (A.55) 



We will show that 

Indeed, assume without loss of generality that i ^ k,l. Then by (|3,52p fij is odd in the 
variables z m i,m € Ai defined in (|3.19p . One the other hand the function &a ,m Rkifki&a' m is 
even in these variables (Pki, HA k , Ha u fkl ac t on different variables). Hence (IA.55|) follows. 
Using that (IA~52|) implies that R^fifia'™ = Ylk^j <t> A ?' ' R ijfi^AT ^Api we have > for 
ij = kl, 

M = <MM +*L > R - f f A * Al ■ (A.56) 
The last relation together with (|A.24|> . ()A.50|) and (|A.55|) . implies the relation (|A.26p . □ 
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B Beyond Born-Oppenheimer Approximation 

Let ipBo(x,y) be the ground state of H^iy), normalized as f \tpBo(x, y)\ 2 dx = 1. Using the 
Feshbach - Schur map with the orthogonal projection 



(Pf)(x,y) = ifcio(x,y) I ipBo{x,y)f{x,y)dx 

(i.e. integrating out the electronic degrees of freedom), as it is done e.g. in [GSJ, we can show 
that 



A G a d (H mo i) <H- A € a d (H nuc i(X)), 

with the corresponding eigenfunctions related accordingly. Here H nuc \(X) is the operator on 
L 2 {R 3M ), defined by 



^nud(A) --=-J2 ^- A vi + E *(y> A )' 

j= i m o 

where k := m/ min^ rrij and 

M 

£ K (y,A) =£(y) + ^— - / Wy^Bo\ 2 dx + 0(k 2 ), 
i zm i J 

with 0(k 2 ) standing for a non-local operator of the indicated order. The energy E K (y) (not 
potential anymore) can be used to define the interaction energy in all orders. 
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